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TECHNICAL NOTE 3350 


THE LINEARIZED EQUATIONS OF MOTION UNDERLYING THE DYNAMIC 
STABILITY OF AIRCRAFT , SPINNING PROJECTILES, 

AND SYMMETRICAL MISSILES 
By A. C. Charters 


SUMMARY 


Linearized equations of motion are derived for both conventlonal 
aircraft having mirror symmetry and spinning projectiles or missiles 
having rotational and mirror symmetry. The aerodynamic coefficients are 
introduced as a formal series expansion in the customary variables, with 
additional terms being included to account for the aerodynamic effects 
of spin. The requirements of symmetry are used to reduce the system of 
aerodynamic coefficients and, in the projectile or missile case, to clarify 
markedly the geometry. A common mathematical approach and standard NACA 
nomenclature are used throughout. 

The equations for aircraft are compared with those for missiles and 
shortcomings in the currently accepted theories are pointed out. The 
dynamic-stability requirements for spin-stabilized projectiles are dis- 
cussed briefly. 

The results are applied to the analysis of flight-test data from the 
aerodynamics range. Relations are derived between the aerodynamic coeffi- 
cients and the constants of the equations of motion. A comparison is made 
with ballistic theory in current use and is found to be satisfactory. 


INTRODUCTION 


The dynamic stability of aircraft is a subject that has been explored 
at great length. The theory of the motion of spinning projectiles has 
also been studied extensively. No further development of either subject 
is needed per se, at least in so far as first-order effects are concerned. 
On the other hand, modem trends in aeronautics and ballistics have brou^t 
the flight performance and even the physical appearance of aircraft and 
projectiles closer and closer together. Hence, a need now exists for a 
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theory of motion that covers both cases simultaneously ivlth a single nomen- 
clature and a common mathematical development. 

In the present analysis^ aircraft are distinguished from projectiles 
and symmetrical missiles as follows: 

1. Aircraft are assumed (a) to have a plane of mirror symmetry 
through the longitudinal axis and (b) to fly only slightly disturbed from 
a steady-state equilibrium attitude so that all components of the angular 
velocity of the aircraft are small. 

2. Projectiles and symmetrical missiles are assumed (a) to have not 
only a plane of mirror symmetry but also 90 ° rotational symmetry (or its 
equivalent) and (b) to fly similarly to aircraft except that the axial 
component of the angular velocity the spin^ may be large (with the 
restriction that the change in the spin must be small). 

Despite these differences^ the flight of both aircraft and projec- 
tiles takes place under closely similar circumstances. The analysis of 
the motion in each case is the classical analysis of a rigid body moving 
under the action of external forces. The differential equations of motion 
are derived in both cases from the vector equations relating the rates of 
change of the linear and angular momenta to the external force and moment. 
Furthermore^ the conditions of f3J.ght prescribed are the -same^ namely^ 
both must fly in nearly a straight line at nearly constant velocity and 
the inclination to the flight path must be small. The aerodynamic-force 
system postulated is the same: The aerodynamic forces and moments are 

assumed to be linear functions of the velocity and the angular velocity. 

In both theories the equations of motion are linearized by the neglect of 
second-order terms. One might well suppose that the dynamic stability of 
aircraft and projectiles had been treated by a common development. Such 
is far from the case. 

The dynamic stability of aircraft was first analyzed by Lanchester 
around I 9 OO; the corresponding analysis of the dynamic stability of spin- 
stabilized projectiles was made by Fowler and his' associates in 1920 (see 
refs. 1 and 2). Both men had very practical objectives in mind. Lanchester 
was interested in the flying qualities of airplanes. Fowler wished to 
find the design criteria for an artillery shell that would insure a true 
flight to the target and a strike head on. Both attacked the problem inde- 
pendently and their divergence of interest is strongly reflected in the 
two developments of the theory^ a divergence which has continued for all 
practical purposes to_ the present time. 

A casual observer reading the literature on dynamic stability would 
be left with the Impression that ballisticians and aeronautical engineers 
were concerned with entirely different problems. The nomenclature is 
different; the geometry does not appear at first glance to be related; the 
mathematical treatments differ radically. The casual observer’s judgment 
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would not be so superficial after all, for it is only fair to say that 
the theory developed to -handle the dynamic stability of aircraft is not 
adequate in its present form to predict the dynamic stability of spinning 
projectiles, and vice versa. For example, the aircraft equations do not 
describe the gyroscopic nutation and precession of a spinning shell, and 
the projectile equations do not describe the phugoid oscillation of air- 
craft. 


Recently, R. E. Bolz and J. D. Nicolaides have derived the dynamic 
stability of spinning projectiles and symmetrical missiles in terms famil- 
iar to the aerodynamicist (see refs. 3 and 4). Although their derivations 
go a long ways toward joining the theories of Lanchester and Fowler, both 
omit the force of gravity and contain certain other restrictions with the 
result that a gap as yet remains between the case of the aircraft, flying 
with its weight balanced by lift, and the case of the gyros copically 
stabilized projectile, flying with varying velocity and spin. It is the 
purpose of this paper to bridge the gap, that is, to treat the two cases 
with strictly similar mathematical developments and with a common nomen- 
clatirre. In addition, it is desired to apply the results of this theory 
to the analysis of experimental data obtained from a relatively new flight- 
test facility, the aerodynamics range. 

The development of the theory presented herein follows closely the 
customary treatment of the dynamic stability of aircraft with the primary 
difference being the simultaneous treatment of longitudinal and lateral 
stability. The nomenclature conforms throughout to NACA standards. The 
other departures of importance from conventional aircraft theory are a 
rather ■ formal development of the aerodynamic force system and the thorough 
use of the condition of symmetry assumed for the body, as is done in the 
ballistic theory. In fact, it might be said that the present development 
is a welding together of these aspects of the conventional aerodynamic and 
ballistic theories which, in the author's opinion, represent the most 
effective means of attacking the problem. 


SYMBOIB 


S'S J 


constants in the x(t) equation 


A 


Ai, Aa 


any quantity 

(This symbol is used in the development of general transforma- 
tions. ) 

constants in the equation 


^o^ 



constants in the 


A(x''') equation 
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^Q} ^X> 

C 

C 


Cm 

<^4 

D 


ID) 


^ot ^x> 
^ 2 > ^3 

F 

g 

H 

i 

1 

I 

Cxz 


K 


constants in the cp(x^) equation 

force 


moment 


aerodynamic coefficients: ; - — — 

(p/2)v2s (p/2)V^SZ 

coefficient of aerodynamic asymmetry force 
coefficient of aerodynamic asymmetry moment 

constants in the |(x'*') equation 
dA 

operator^ DA = -- 

CL Li 


operator j IDA = 


Cpig + CjTg 


dA 


dx‘ 


constants in the t](x+) equation 

force, external, acting on body 
acceleration of gravity 
angular momentum of body 
unit vector along X axis 

moment of inertia of body 
unit vector along Y axis 
product of inertia of body 
unit vector along Z axis 


lY 


Kv 


i2L 

mZ^ 
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Kz^ 


Ki 




Ix. 

mZ^ 

iL 

mZ^ 

£h 

% 

■^XZ 

iz 


z 

b 

'P 

Z. 

P 


ctiaracteriBtic length 


4Kx2 

i^Kx2 

2 Kx^ 





!fl£ 

2 Kx^ 


m 

mo 

mo 

mq 


mass of body 

P-(Cm3 - 1Cm2^ 
2 Kz^ 

4Ky^ 

Cmq_ 
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ma 


”1 


2 Ky^ 

Cm. 

g 

4Ky^ 

mq + Xq + iv(K - mj^-p) 

ma + ivmpp 

^■d ’*' 


“rp 


Cmj<p 

8(iKy^ 

Sp 

4Ky^ 

! 3 e 

8pKy^ 


M 

M 

np 

nr 

np 

n. 

P 


linear momentum of body- 
moment j external acting on body 
^np 

4 Kz^ 

4K/ 

2Kz^ 

kKz^ 
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•'ll. 




nqp 


°ccp 


■‘op 


qp 




SpKzf 


P 

P' 

Po 

P(S) 

q. 

Q(S) 

r 

R(S) 


S 

S 

t 


X component of cj 


variation in p over measured trajectory 


constant component of angular velocity about ■which, the angular 
velocity in flight varies (p = p^ + p’) 

polynomial factor in cons-tants dg of iCx"*") equation 


Y component of 10 


polynomial factor in cons-tants f2, fs of Ti(x^) equation 


Z component of co 


stability quartic 


stability factor. 


w 


arc length along trajectory 
characteristic area 

independent variable in stability quadratic 
roots of stability quartic, Sx = dg, ^2 = d^ 

time 
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Ti, Ts terms In SxR anti SgR 

u variation in velocity along trajectory 

V velocity of center of gravity of body 


constant component of velocity about which the flight velocity 
varies (V = Vq + u) 


3 j + ak 


V 

angular velocity of body with respect to XYZ axes 


X 

distance along Xq axis 


x+ 

X 



pi 



mg sin 7 q 



PSVqS 

7 


^Xo 


Xo 

2 



^Xg 



2 


X. 



a 



Xn 

C-y- 


4 . 






^oD 

2 

*- 





2 

■f 
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m 

Xpp 

Cv 

■'^PP 

8p^ 




X 

coordinate axis 





Xo 

space-fixed coordinate axis 




1 

distance along Yq axis 




it 

y 

pi 





Cy 

2 




Y* 

% 




4p 




7p 

Cy 

Xp 




Xr 

^Yr 




XpL 

Xp + ^0 




Xxp 

^Xccp 

V 




^dp 

Cy* 

■^ap 

8p2 




Xqp 

‘^Xqp 

8p2 




Y 

coordinate axis 




Yo 

space-fixed coordinate axis 




z 

distance along Zq axis 





z 

pi 
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Za 

z. 

a 

"Tl 

ZoL 

ZaL 

^rp 

Z 

Zo 


mg cos 7 q 


1 - Zq + Ivz, 


^ + iVZpp 


1 + z. + ivzi 

a PP 


, 

~2~) also^ ^aL 


Za - Xq 


coordinate axis 


space-fixed coordinate axis 
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angle of attack, 90° - ZV 
angle of sideslip, 90° - YV 


inclination of fligjit path of steady-state glide to 
horizontal, ^ gZo 

y + iZ , nondimens ional transverse displacement 


variable of integration, replaces x"** in integrand 


t - ie 


angle of pitch, angle between intersection of XoYq ^nd XZq 
planes and X 


arbitrary angle of rotation about X 


aircraft density factor. 


P + ia 


density of the air 


'l-i 

s 


time factor. 


angle of roll: 

aircraft: angle between intersection of YZ and XqYq planes 

and Y; 

symmetrical missile: angle between Y and 2 

angle of yaw, angle between Xq and intersection of XqYq 
and XZq planes 



12 


NACA TN 3350 


w 

angular velocity of body with respect to space-fixed axes 

- 

a 

angular velocity of 
axes 

the XYZ axes with respect to the 


? 

1, 2, 3 

body-fixed coordinate axes, fixed in body for missiles 
projectiles 

and 

— 


angle, ^ XY = angle between X and Y 


- 

(*) 

operator, A = — 
dt 




(••) 

operator, A = — & 
dt^ 

— 



n 

a vector quantity 

Subscripts 


*• 

The definitions of all subscripts apply to their usage with all sym- 
bols except where the complete symbol with its subscripts is defined as 
a unit. 


A 

aerodynamic components 



D 

component along trajectory (drag) 



G 

gravity force components 



I 

imaginary component 


- 


1 

component along X 

axis 


- 

L 

component normal to 

trajectory (lift) 



m 

component along Y 

axis 


-- 

n 

component along Z 

axis 


* 

P 

operation, A-n = ^ 
op 



I 

q 

? 1 A 

operation, = •^ 



fa 
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R 

X 

Y 

Z 

a 


operation, A~ = 

or 


real component 


component along X 


component along Y 


component along Z 
operation, A^ = 

oa 


a 


operation, A<v = 

oa 


operation, Aj 3 = 




0 

1 
a 


operation, ^ 

value of coefficient witli u = p 

value of quantity at z = 0 or x 
component along 1 axis 
component along 2 axis 
component along 3 axis 


= a= P= co = p = q = r = 0 


= 0 


EQUATIONS OP MOTION FOR AIRCRAFT WITH 
ONE PIARE OF MIRR(® SYMMETRY 


The essential features of Lanchester’s work have formed the basis of 
all subsequent treatments of this subject, although many details such as 
nomenclature have been modified to meet better the changing demands of air- 
craft design. Dynamic stability is a standard subject in aeronautical 
textbooks (e.g., see refs. 1 and 5 ) and, since the methods and results of 
the theory are well known, the aircraft case will be sketched rapidly. 
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Geometry 


The geometry is shown in figures l(a), l(b), and l(c). Two sets of 
Cartesian coordinate axes are used in developing the equations of motion. 
One set is fixed in the aircraft and is known as "stability axes." The 
other set is fixed with respect to the earth and is known as "earth axes." 
The orientation of the two sets is shown in figure l(a). Figure l(b) 
shows the location of the gravity vector with respect to earth axes and 
figure l(c) the location of the velocity vector with respect to stability 
axes . In the figures ■ the directions of the coordinate axes are indicated 
by unit vectors and the angles by circular arcs, which represent great 
circles on a unit sphere centered at the center of gravity. Both sets 
of axes are right-handed, as defined in reference 1 ; that is, the positive 
sense of a component rotation or couple about any axis is "determined by 
reference to a right-handed screw, when facing the positive direction of 
the axis." It will be noted that the various features of this geometry 
conform throughout to NACA standards. 

The stability axes are designated by XYZ. They are fixed in the 
aircraft and are oriented so that Y is perpendicular to the plane of 
symmetry and Z is perpendicular to the relative wind in steady flight. 
Their origin of coordinates is at the center of gravity. In other words, 
they are body-fixed axes with the XZ plane being the aircraft’s plane 
of symmetry and with the X axis being coincident with the velocity vec- 
tor in the steady-state glide on which the flight path is a perturbation. 
The positive directions along the axes are as follows: X is positive 
forward; Y is positive to starboard; Z is positive towards the bottom 
of the aircraft. It should be noted that, whereas Y is always a prin- 
cipal inertia axis, X and Z are not necessarily sc, and products of 
inertia terms are Introduced into the equations of motion in order to 
accoiuit for the alinement of X with the steady-state relative wind 
rather than with a principal inertial direction (body axes alined in this 
particular .way are known as stability axes). 

The .earth axes .pre designated by XqYoZo. They are fixed with 
respect to the earth, but, since the earth’s rotation is negligibly slow 
compared to the angular velocity of the aircraft, they are considered to 
be Galilean axes fixed in space. As with the stability axes, their ori- 
entation is determined by the position of the aircraft in. its steady- 
state glide, which is subject to the restriction that in the steady-state 
glide the plane of symmetry is a vertical plane. Each member of the 
XqYqZq axes is taken to be parallel to and to point in the same direc- 
tion as its corresponding member of the XYZ axes. In other words, X^ 
points in the direction of the steady-state glide path, Yq is horizontal 
and points to the right (of an observer facing forward), and Zq points 
down (but not necessarily vertically down). Two locations are designated 
for the origin of coordinates of the earth axes, depending on the com- 
ponent of the motion being considered. For angular measirrements , the 
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origin of coordinates is taken to be the center of gravity of the air- 
craft at some particular point along its flight path. For linear measure- 
ments, the origin of coordinates is located as required by the over-all 
scheme of measurement vhich would be set up to analyze the flight of the 
aircraft. 

The orientation of the aircraft in space is defined by three angular 
coordinates cp, 0 , and ijr which give the alinement of the stability axes, 
XYZ, with respect to the earth axes, XqYoZo (see fig. 1(a)). If these 
angles are all zero, the stability axes point in the directions of the 
earth axes. Any other orientation of the aircraft is reached by three 
consecutive rotations, starting with XYZ pointing in the directions, of 
XqYqZq. In defining each of the rotations, the point of view is taken 
that XYZ are either in their starting position or in the position given 
by the preceding rotation, and not necessarily in their final position. 

The rotations are noncommutative and must be taken in the order specified. 

ijr Rotation 1. Start with XYZ pointing in the directions of 

XqYqZq, rotate about Z through Ic. 

6 Rotation 2. Rotate about, Y through 0, thereby bringing 

X to its final position. 

<p "Rotation 3. Rotate about X through cp, thereby bringing 

Y and Z to final positions and the aircraft to its actual 
orientation in space. 

The angles "'1^, 0, and 9 may also be defined as the angles between 
coordinate axes and the intersections of certain planes. For the purpose 
of these definltiona it is considered that the origin of coordinates of 
the XqYqZq axes is momentarily coincident with that of the XYZ axes . 
The point of view taken here is that the XYZ axes are in their final 
positions given by the orientation of the aircraft in space at the moment 
in question. The sign of the angle is specified by giving the axis about 
which the rotation is taken in going from the line named first to the 
line named second in the definition. The axis of rotation is listed in 
parentheses in each definition. 

angle between Xq and intersection of XqYq and XZq planes 
(rotation about Zq) 

0 angle between intersection of XqYq and XZq planes and X 
(rotation about Y^' = Yq rotated about Zq through t^) 

9 angle between intersection of YZ and XqYq planes and Yj also 
angle between Intersection of YZ and XZq planes and Z 
(rotation about X) 
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The inclination of the earth axes (XqYqZq) to the horizontal is 
given by the angle 7 q (see fig. 1 (b)). The angle 7 q is positive when 
Xq points above the horizontal (as shown in fig. 1 (b)). It should be 
noted that Yq is horizontal by definition and hence the XqZq plane is 
a vertical plane and contains the gravity vector 'g. 

The orientation of the velocity vector V with respect to the sta- 
bility axes XYZ is given by the angle of attack a and the angle of 
sideslip p (see fig. 1 (c)). These angles are defined so that the product 

of the velocity magnitude V and the sine of p or a gives the component 

of V along Y or Z; that is^ by definition 

p = 90° - 

a = 90 ° - ^YZ 

and resolving V along the transverse stability axes YZ results in the 
desired relations, as follows 

Vy = V cos VY = V sin P 

V2 = V cos ^ VZ = V sin a 

The above definitions of a and p were "chosen because it is believed 
that the transverse components of the velocity are the quantities having 
the most physical significance insofar as the aerodynamics of the air- 
craft is concerned. However, it should be noted that these definitions 
differ somewhat from common urtnd -tunnel practice. In a wind tunnel the 
model is normally placed at an angle of sideslip by supporting it on a 
bent sting. The angle of attack is set at zero, and the sting is alined 
to lie in the plane containing the axis of rotation of the angle-of-attack 
sector and the tunnel axis (it is assumed that the wind vector is along 
the tunnel axis). The model is then placed on the sting with its plane of 
symmetry perpendicular to this plane. The angle of attack is now varied 
by rotating the angle-of-attack sector. Consequently, the Z axis of the 
model rotates in a plane perpendicular to this sector's axis of rotation. 
The angles of attack, a, and sideslip, p, are defined as follows: a is 

the angle through which the angle-of-attack sector is rotated; p is the 
angle at which the sting is bent. Now, referring to figure l(c), it is 
clear that a defined in thl£ report is the same as a defined in wind- 
tunnel practice, since the ZY plane shown in this figure is the plane 
in which Z rotates and a is measured in this plane i^n both cases. On 
the other hand, p in this report is measured in the YY plane whereas 
P in the wind tunnel is measured in the XY plane, as shown. The dif- 
ference is second-order for small values of p and hence is insignificant 
in so far as linearized theory is concerned. However, should the theory 
be extended to the case of aircraft maneuvering at large angles, as is 
done in reference 6 , the precise definitions of a and p become important. 
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It will be convenient to write down expressions for the direction 
cosines between various axes of the earth axes and stability axes systems. 
The direction cosines are readily obtained from the equations derived in 
section (5.0)^ "Kinematics/' of reference 6 and are listed in the table 
below; each entry in the table gives the direction cosine between the axis 
heading the row and the axis heading the column belonging to the space in 
which the entry is listed. 



X 

Y 

Z 

Xo 

cos 6 cos 4" 

cos sin T sin 0 - 

sin lie cos 

cos ijf cos cp sin Q + 
sin ijf sin (p 

Yo 

cos 9 sin ijf 

sin sin <p sin 6 + 

cos ijf cos q> 

sin ijf cos cp sin ,0 - 
cos ijf sin cp 

Zo 

- sin 0 

sin cp cos 6 

cos 9 cos cp 


Now, the conditions postvilated for the flight, which are discussed in the 
next section, require that conresponding axes of the two coordinate systems 
lie close to one another during the flight covered by the analysis and 
hence that the angles cp, 0, and ''I'' be small. Only first-order terms are 
retained in the development of the equations of motion and, since cp, 6, 
and are first-order, quadratic terms in the expressions for the direc- 
tion cosines will be second-order and may be neglected. The equations 
giving the direction cosines coirect to first order are summarized in the 
table below; again each entry in the table gives the direction cosine 
between the axis heading the row and the axis heading the column belonging 
to the space. in which the entry is listed. 



X 

Y 

Z 

Xo 

1 


e 

Yo 


1 

-9 

Zo 

-e 

9 

1 


( 1 ) 
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Conditions of Flight 


The folloving conditions are postulated for the aircraft and for its 
motion in flight: 


(1) The aircraft has a plane of mirror symmetry, the XZ plane. 

(2) The magnitude of the velocity of the aircraft may be large but 
the change in the velocity must be small. Accordingly, the velocity may 
be represented by a constant plus a small perturbation; that is, V=Vo+u. 

(3) The angles orienting the two sets of coordinate axes, Q, and 
and locating the velocity vector, p and a, are all small. 

(4) All components of the angular velocity, p, q, and r, are small. 

(5) The flight path of the aircraft is a perturbation on a steady- 
state glide. 

Conditions (2), (3)^ and (4) may be summarized by saying that the 
following variables must be small (e.g., have a numerical magnitude of 
0.1 or less) and hence are first-order quantities: 


9 e h 


-iL p 
Vn 


a 


qZ rZ 
2Vo 2Vo 2Vo 


small, for example, 
0.1 or less 


( 2 ) 




where Z is a characteristic length. ^ It should be noted that if P; 
2^ ^Z PZ_ d'Z 


order, where B = -7^ and a 
^ dt 


are first-order, then 
da 
dt * 




2Nr 


will also be flrst- 


The theory developed in this paper includes only first-order terms . 
Products of first-order terms ( except for the Magnus terms of spinning 
projectiles and symmetrical missiles), for example, 90, are considered to 
be negligibly small and are left out of the equations of motion. As a 
result, the equations of motion become linear in form and their solution 
may be given by well-known, explicit functions. Experience has sho™ 
that the solution of the linearized equations of motion is, in fact, a 
reasonable representation of the flight motion of many aircraft. However, 
there may be circumstances under which the linearized equations do not 
describe the aircraft's flight, despite the fact that all the quantities 
listed in equation (2) are small. The aerodynamic characteristics of 
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certain aircraft are such that some of the quadratic terms are as large 
or larger than certain of the first-order terms included in the linear 
theory. For aircraft of this type, the significant quadratic terms must 
be retained in the equations of motion and the theory revised accordingly. 
Consequently, the analysis of an actual flight should include a check of 
the magnitude of second-order terms to be sure that they are in fact 
second-order, that is, negligible. In other words, the a priori assump- 
tions of the theory should be checked a posteriori by an analysis of the 
measurements . 


Kinematic Relations 


The equations of motion are derived from the basic vector equations 



( 3 ) 

dt 

^ = M 
dt 



The first step in their derivation is to obtain relations between 
the physical quantities, M and H, and the kinematic variables (of the 
motion), 3 , a, cp, 9 , and . The components of the linear momentum, M 
along the XYZ axes in terms of the components of the velocity vector, 
V, and the components of V in terms of p and a are listed in the 
following table. 



M 

V 

i 

mVx 

Vx = V^/■l - sin^a - sin^P 


mVy 

Vy = V sin P 

k 

mVg 

V2 = V sin a 


( 5 ) 


where the component along X, Y, or Z is designated by i, j, or k, the 
unit vector along X,Y,Z, respectively. 

Similarly, the components of the angular momentum, H, are given in 
terms of the components of the angular velocity, m, and, in turn, the 
components of to in terms of given in the table below 

(see ref. 6 ) . 
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H 

0) 

i 

- '^XZ^ 

• • 

p = cp - ijr sin 9 

J 

lyq 

a • 

q = 0 cos cp + t sin cp cos 6 

k 

l2r - Jx^ 

r = li' cos 9 cos cp - 6 sin cp 


( 6 ) 


7 


For convenience In computing the relative magnitudes of terms in the 
equations of motion, the components of w (p, q, and r) may be written 
as follows, expanding sines and cosines and retaining only first- and 
second-order terms. 


p = (p - ar 1 


q = a + cpr 


(7) 


r = - cpq j 

The expressions for V and w given in equations ( 5 ) and ( 7 ) include 
second-order terms in order to clarify certain steps taken in linearizing 
the equations of motion. The components of V and io appearing in the 
inertial terms in the equations of motion are multiplied by large numeri- 
cal factors. Consequently, it is desirable to compare the second-order 
inertial terms with the first-order aerodynamic terms in order to estab- 
lish that the second-order Inertial terms are truly negligible. The com- 
parison is made in Appendix A, which covers the linearization procedure. 


Physical quantities are defined with respect to the body-fixed axes, 
XYZ, as is customary in the treatment of the motion of a rigid body. On 
the other hand, the rates of change Involved in the basic vector equa- 
tions ( 3 ) and (4) are taken with respect to space-fixed axes. Hence, 
equations ( 3 ) and (4) must be transformed in order to account for the 
movement of the XYZ axes . The well-known transformation for the rate 
of change of any vector A from fixed to moving axes is given by 


— = A + w X A 
dt 


( 8 ) 


where A is the rate of change of A measured with respect to XYZ 
axes. The components of the transformation are 

d^ 

dt 


dAy * 

= Ax + qAz - rAy ^ 


dAy 

— = Ay + rAx - V^z 
^ = A 2 + pAy - qAx 


( 9 ) 
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Aerodynamic and Gravitational Force Systems 


The resultant external force, F, of equation (3) Is the sum of the 
aerodynamic forces and the gravitational force; that is, 

F = F^ + Fq ( 10 ) 

The resultant external moment, M, of equation ( 4 ) is due solely to aero- 
dynamic reactions . 


Concerning the aerodynamic forces and moments, it is assumed that the 
components of the resultant force and moment, F^ and M, are given by 


^Ax,y,z = ^x,Y,z I 


( 11 ) 


and 


%,Y,Z “ ^Z,m,n 2 




(12) 


It is further assumed that the coefficients, the C*s, in equations (ll) 
and (12) are functions of the variables cc,, p, <x, p, q, and r and that 
these functions may be expanded as a series in the variables named. The 
general formula for any coefficient, Cq^, where a stands for X, Y, Z, I, 
m, or n, is 


Co = C, 


+ C; 


apP + C, 




pi 

,0- + Ca_. ^ ^ 


r 0,1 


^ + ^^q 2V ^ 


r + 

% 2 V ^ 


Caj, + (higher-order terms) (l 3 ) 

The coefficients of this series are assumed to be independent of the 
variables P, a, A, p, q, and r and to be functions only of the exter- 
nal shape of the aircraft and of fundamental aerodynamic parameters , such 
as the Reynolds and Mach numbers . They will be known herein as the 
"aerodynamic coefficients" and correspond to the conventional stability 
derivatives . 

Examples of the higher-order terms are 
quadratic terms of the series. In the aircraft case all of the variables 
P, a, p, q, and r are small (first-order). Hence, their products will 
be second-order and the quadratic and higher-order terms of the series may 
be neglected. 

If the aircraft had no plane of symmetry, all 48 coefficients would 
be required to account for the complete force and moment system. However, 
the majority of aircraft have one plane of mirror symmetry as postulated 
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i 

here and^ aa a result, certain of the coefficients must be zero. The 
requirement imposed by the existence of a plane of mirror symmetry can be 
stated succinctly as follows. 

Let the aircraft execute two motions, the first being the mirror 
image of the second about the XZ plane, the plane of symmetry; the aero- 
dynamic force due to the first motion must be the mirror image of that 
due to the second. The mirror image of a motion (or force) is defined as 
the motion (or the force) that will look the same to a second observer 
stationed on the opposite side of the aircraft from the first; that is, 
one observer being to port, the other to starboard. 

To illustrate, consider the terms C2p(pZ/2V), the rolling moment due 
to rolling velocity, and Cmp(pZ/2V) , the pitching moment due to rolling 
velocity. The mirror image of the rolling velocity p is always -p. 

The mirror image of the rolling moment is -Mjr; hence, the mirror 

image of C2p(pZ/2V) should be -C2p(pZ/2V) . It can be seen that this term 
reverses sign as it should, since p mirrors to -p, and the coefficient 
Cl is therefore ah allowable one. On the other hand, the mirror image 

Jr 

of the pitching moment My is -riMY; hence, the mirror image of C]]jp(pZ/2V) 
should be +Cinp ( p Z /2V ) . It can be seen that the term reverses sign as it 
should not, since p mirrors to ^p, and therefore the coefficient 
is not an allowable one. Consequently the coefficient Cm_ must vanish; 
that is, Cmp =0. ^ 

Applying the criterion of symmetry to equation (l3), one’ can see that 
the following coefficients must vanish: 


Cy 

^3 

= Cy. 

= Cy = 
Xp 

0 

>c 

II 

0 '' 

0 

h; 

0 

= Cy 
■“•a 

- Cy. = 
■••a 

% “ 

0 



= = 


0 


= ^Z 

= = 
a 

Cj * 

0 

Cmp 

P 

*"mp = 

Cmr - 

0 




^nq = 

0 


> 


(14) 


Substituting equations (lU) in eqiiation (l3) and arbitrarily assigning 
a minus sign to the coefficients of Cy and in conformity with standard 
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practice (since the axial drag force usually acts in the -X direction 
and the normal force in the -Z direction), give the force and moment 
coefficients for the aircraft as follows; 


Cx 


aZ 


-CXo - ^ " ^X, 




•Q, 2 V ^<1 2 V 




Cy = + Cyj ^ Cy — — + Cy 

i ^ ^ 2 V 2 V 2 V 

Cz = -CZo - Cz^a - ^ 

= ^ZpP + ^Z^ 2V'*'‘^^r ^ 


(15) 


*^n ^n^ ^ + *^np 2 V "** 2 V 


Concerning the force of gravity, Fq, Its components along the XYZ 
axes are readily computed by first resolving Fq along the XqYqZq axes 
and then computing the XYZ components using the direction cosines given 
in equations (l). The force of gravity Fq lies in the X^Z^ plane 
since Yq is horizontal by definition and, hence, is given by 


Fq = lo(-mg sin 7 q) + ko(mg cos 7^) 


(16) 


The XYZ 
components of 


components of Fq are now obtained by resolving the XqZq 
equation (16) along the XYZ axes using equations (l); 


N 

= mg (- sin 7o - 0 COB 7 q) 

Fq^ = mg (ijf sin 7 o + 9 cos 7 q) 

Fq = mg (- 0 sin 7 q + cos 7^) 


(IT) 


Derivation of the Differential Equations of Motion 


The differential equations of motion are derived from the basic vec- 
tor equations (3) and (T) by resolving the vectors along coordinates axes 
and using the kinematic relations and the equations for the external 
forces and moments, thereby obtaining six scalar equations, three asso- 
ciated with the force and three with the moment. Both earth axes and 
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stability axes are employed as coordinate axes in resolving tlie basic 
vector equations. Earth axes are used to compute the motion of the center 
of gravity of the aircraft, since linear displacements are measured along 
these axes. Stability axes are used to compute the angular motions of 
the aircraft and of the velocity vector, since the inertial properties and 
aerodynamic coefficients of the aircraft are invariant with respect to 
these axes. The detailed steps taken in deriving the equations of motion 
are outlined as follows . 

Both equations (3) and ( 4 ) are transformed from fixed to moving axes 
by equation ( 9 )^ except for the X component of equation ( 3 )^ as dis- 
cussed below. Equation (l?) for* the gravity force and equations (U) and 
(15) for the aerodynamic forces are substituted into equation (lO) for the 
external force. The kinematic relations for linear momentum and velocity 
are given by equations (5) and these are used to complete the development 
of equation (3). Equations (l 2 ) and (15) for the aerodynamic moments and 
equations ( 6 ) and (?) giving the kinematic relations for the angular momen- 
tum and angular velocity are substituted into equation ( 4 ) to complete its 
development. 

The details of this particular step are given in Appendix A. It 
should be noted that the equation associated with the axial drag force is 
treated as a special case. This eq-uation is derived by resolving the 
vector linear momentum equation ( 3 ) along the space-fixed Xq axis rather 
than the X axis, thereby avoiding the term (qA2 - rAy) which comes from 
the transformation, from fixed to moving axes (eq. 9 )« Th® change in 
approach is desirable because in the linearization of the equation of 
motion the term (qA^; - rAy) gives rise to second-order terms whose magni- 
tude compared to first-order terms is imcertain and whose neglect conse- 
quently may be questionable. In so far as the validity of the derivation 
is concerned, either axis system, stability axes or earth axes, may be 
used in deriving any of the equations of motion provided the vector quan- 
tities are resolved correctly along the axes in question and the angular 
velocity of the coordinate axes (if any) is accounted for. 

The next step in the derivation is to separate each of the six epila- 
tions into first- and second-order terms (and higher-order terms) and to 
neglect the second-order (and higher-order) terms. This step is described 
in detail in Appendix A. Hence, at this point, the equations of motion . 
are linearized and the theory is restricted to a consideration of first- 
order effects. 

The third step in the derivation is to introduce a nondlmensional 
time as the independent variable and a concise notation for the factors 
of the dependent variables. This nondlmensional time, f*", is sometimes 
referred to as the reduced time and is based on the aircraft time factor, 

T. The concise notation stems directly from the subscripts of the aero- 
dynamic coefficients and at the same time includes any other multiplica- 
tive factors. For example. 
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% 


2Ky^ 


In the same spirit^ the derivative notation, is shortened by defining 


DA 


dA 

dt+ 


This step parallels a similar procedure followed in reference 7 * The 
advantages accrued thereby speak for themselves. 


The final step goes back to (5) of the condition postulated for the 
flighty namely^ that the flight path is a perturbation on a steady-state 
glide. Hence^ the equillbriiun condition of flight with all the dependent 
variables (ujP,a^p,o,^p,qjr) equal to zero is one in which the lift, drag, 
speed, glide angle, and weight are all in balance. It is assumed here 
that the Xq axis and the X axis are both coincident with the trajec- 
tory along the steady-state glide path. The equations of motion for the 
steady-state glide are 




mg sin 7 q 

pSVo^ 


z 


o 


mg- cos 7 q _ 
pSVo^ " 


> 


(18) 


mo = 0 


Equations (I8) are used to replace the gravity terms -with corresponding 
aerodynamic terms in the final differential equations of motion. 


At the same time, the factors replaced by the 

factors XqJ),XqJ)^ 2 ol,Z(^P which refer to the drag and lift of the aircraft 
as measured in a wind -tunnel. Now, it is cus-tomary in deriving the equa- 
tions of motion to use drag and lift coefficients rather than the X and 
Z coefficients. This is a perfectly reasonable choice since the static 
s-bability characteristics of many aircraft are measured in wind--tunnel 
tests. Wind--tunnel data are ordinarily presented as lift and drag polars 
and it is desirable to formulate the theory in terms of lift and drag so 
that these data can be used directly. A support system commonly employed 
in wind -tunnels has already been described. In so far as the coordinate 
axes of this report are concerned, lift and drag are measured along and- 
perpendicular to the steady-s-bate relative -wind, that is, along the -Zq 
and -Xq axes, respectively. The aerodynamic lift and drag factors 
corresponding to the lift and drag coefficients are designated by the 
subscripts L and D, and the relations between these and the aerodynamic 
factors corresponding to the X and Z coefficients are given by 
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= ^oD 




- 2qL 


= Z T 

o oL 


^cc = + ^oD 


( 19 ) 


The six differential equations of motion which result from the pre- 
ceding steps are given below in the following order: the three force 

equations, axial drag, side force, and normal lift force, and the three 
moment equations, roll, pitch, and yaw. 


^VVo/ 

) + 2x^jj 

( u \ 

x^Da + 

- 

Zol)cl + x;gD0 

+ Zop0 = 0 

( 20 ) 

(1 - 

y-) Dp - 

p 

73^ + 

(1 - 7 r 

) Dt + 

"^oD^ - yp°^ - 

Zol9 = 0 

(21) 

(1 + z^) 

Da + ( 

+ ^oD^ 

)a - (1 

- Zg) 

D0 - x^p0 ,+ 2z 


(22) 


D^9 - 

2pDCp - 

KiD^^if - 


- Z^DP - ZpP = 

0 

( 23 ) 



D^0 . 

- mqD0 - 

m»Da 

w 

- mQ_a = 0 


(24) 


D^t - 

npDt - 

KgD^cp - 

npDcp 

- n^Dp - npP = 

0 

( 25 ) 


Discussion of the Differential Equations of Motion 


An inspection of equations (20) through ( 25 ) shows that the six equa- 
tions may be divided into two distinct groups of three equations each. 

One group consisting of equations (20), (22), and (2^) involves only the 
variables u/Vq, a, and 0 and describes the longitudinal motion of the 
aircraft. The other group consisting of equations (21 ), (23)^ and ( 25 ) 
Involves only the variables 3, cp, and ilr and describes the lateral motion 
of the aircraft. Hence, it may be concluded that the linearized equations 
of motion of aircraft may be divided into two separate parts and that the 
parts do not interact. This result is, of course, well known and the 
over-all stability of aircraft is customarily subdivided into longitudinal 
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stability involving u/Vq, a, and 0 and into lateral stability involving 
cp^ and i|f. In this connection, it should be noted that different rep- 
resentative lengths, I, may be used in the two groups of equations since 
the groups are independent of one another. In practice, it is common to 
use the mean aerodynamic chord of the wing (z = c) in the longitudinal 
stability equations and the wing span (z = b) in the lateral stability 
equations. 

Equations (20) through ( 25 ) agree satisfactorily with corresponding 
equations presented in standard aeronautical texts (see, e.g., eqs, (lO- 89 a, 
b, c) and (ll-3^a, b, c) of ref. 5)* fhe equations here do contain certain 
extra terms not ordinarily found elsewhere but it is believed that all of 
the terms in these equations are required for a complete and consistent 
first-order aerodynamic force system and to account accurately for gravity 
forces. On the other hand, the additions are so small for representative 
aircraft that they do not affect the comparison in so far as practical 
application of the equations is concerned. 

To summarize briefly, it may be stated that equations (20) through 
( 25 ) are the standard equations underlying the longitudinal and lateral 
stability of aircraft. Hence, it has been demonstrated that the approach 
uned in the present development leads to the commonly accepted result. 

The solution and application of the differential stability equations have 
been treated at length and reference is made to the extensive aeronautical 
literature for a discussion of this aspect of the subject (e.g., see 
refs. 1 , 5 , 6 , 'J , 8 , and 9 ). 


EQUATIONS OF MOTION FOR PROJECTIIIS AND MISSILES 
WITH 90 ° ROTATIONAL AND MIRROR SIMMETRY 


In this paper the equations of motion describing the flights of spin- 
ning projectiles and symmetrical missiles are treated simultaneously, and 
the results may be applied to either case. Both projectiles and missiles 
are mentioned in the title because the means of stabilization customarily 
employed in the two cases depends on different physical principles, the 
majority of artillery shell being spin stabilized and missiles being fin 
stabilized. Despite this difference, the motions of both stem from a 
common set of equations and each represents an application of the general 
equations to a particular condition of flight. The terms, projectiles 
and missiles, will be used interchangeably for the sake of succinctness, 
but, although only one may be named, the other will always be implied. 

Projectiles and missiles are considered to differ from aircraft in 
two respects; 

( 1 ) They are postulated to have 90° roll symmetry or its equivalent 
in addition to a plane of mirror symmetry. 
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(2) Ttie component of their angular velocity along the axis of sym- 
metry is alloved to be large] that is, they may roll continuously (or 
spin) around the axis of symmetry. 

It is immediately evident that the freedom of the missile to spin 
Vill lead to difficulties if axes fixed In the missile are specified in 
this case as they were for aircraft. The angle of roll, cp, ^flll become 
large and the components of the gravity force will vary with the sine and 
cosine of cp. Both of these consequences will violate the requirement 
that the dependent variables remain small during the flight, an essential 
feature of a linearized theory. It is clear that the geometry postulated 
in the aircraft case must be modified before it is suitable for the devel- 
opment of linearized equations of motion for projectiles and missiles. 

The 90° rotational symmetry, the addition to the requirement of mir- 
ror symmetry, provides the key to the problem. It will be shown that the 
aerodynamic force system is modified for this case so that the resultant 
force and moment vectors are independent of the roll orientation (to first 
order). Consequently, only one coordinate axis need be fixed in the 
missile, that is, along the axis of rotational symmetry, the longitudinal 
axis of the missile. One of the transverse coordinate axes will be pre- 
scribed to lie in a certain space-fixed plane oriented so that the angle 
of roll of the coordinate axis system is always zero. As a result, the 
change in orientation of the coordinate axes remains small during flight, 
and the variation of all of the dependent variables except the roll angle 
of the missile is correspondingly small so that the basic requirement of 
a first-order theory can again be satisfied. 

Consequences of 90 ° Rotational Symmetry 


Rinety-degree rotational symmetry may be readily visualized by pic- 
turing the missile (or projectile) in two positions, one rotated by 90° 
about the axis of symmetry mth respect to the other. If the missile has 
90° rotational symmetry, the two pictures will look precisely the same. 

In mathematical terms, rotation through 90 ° transforms the missile into 
itself. It is assumed in this paper that the axis of rotational symmetry 
coincides with the . longitudinal axis of the missile. 

First, the consequences of rotational symmetry to the aerodynamic 
force system will be investigated. Let the XYZ axes remain axes fixed 
in the missile, but place the X axis along the axis of rotational sym- 
metry. The X axis is now a principal axis of Inertia so that 
vanishes, furthermore ly = I2. It is clear that the aerodynamic coeffi- 
cients must have values such that a rotation of the V, Vy, and w vec- 
tors by 90° about X will produce a simiilar rotation of the ^ 

vectors of 90° about X, since F and M are fimctions of V, Vy, and CJ 
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and the aerodynamic configuration will be precisely the same after the 
90° rotation, as required by symmetry. The components of V,j, are given 
to first order, by = Pj + dk. It will be noted that transforms 

as V does. _Only V will be written down explicitly in the subsequent 
analysis but Vij will be included implicitly in all derivations in which 
V, and u are involved. 

Now, if A is any vector and is A rotated 90° around the axis 

of symmetry, and if_the direction of rotation is taken to be from Y to Z, 
the components of A are related to the components of as follows: 



Applying equation ( 26 ) to the 
aerodynamic force and moment vec- 
tors gives the relations between 
the components of the force and 

moment before and after rotation 

• 

of V, Vij, and u, which are 



Sketch (a) 


(V*, = Px (V, 

Fy* (V*. “*) = -Fz (V, oJ) 

Fz* (V*, w,) = Fy (V, aJ) 
(V*, %) = % (V, m) 
My* (V*, = -M2 (V, w) 

Mz* (V*, to^) = My (V, u) 



( 27 ) 


It should be noted that the values of the components after rotation must 
be computed using the rotated values of velocity and angular velocity, 
and as indicated by the fimctional notation of equations (27). 
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Equations (27) must h£ld for all values of V and W. Hence^ com- 
puting the components of and from equations (27) and substituting 

the values of the components of V and w from equations (5), ( 6 ), and 
(7) (neglecting second -order terms) and of and M from equations (ll)^ 
(12)^ and (15)^ the following equalities must hold between the aerody- 
namic coefficients: 



Consequently, the equations for the components of F and M become 
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Second, since the aero- 
dynamic coefficients for rota- 
tionally symmetric bodies have 
been established, the varia- 
tion of the aerodynamic coef- 
ficients with rotation of the 
YZ axes aroxmd X will be 
investigated. Let the YZ 
axes be turned through an 
arbitrary angle A to a new 
position Y^ZT*", as shown in 
the sketch. 

Relations between the 
components of a vector A 
along YZ and along Y^Z^ 
are given by 



Sketch (b) 


Ay = Ay+ cos A - sin A 


= Arr+ 


(30a) 


cos A + Ay'^ sin A 


Ay"*" = Ay cos A + A^ sin 


= A2 cos A 


Ay sin A 


(30b) 


Since the rotation is about X, the X component of A is unchanged. 


Now, vector quantities themselves do not change with rotation of 
coordinate axes; only their components along the axes vary. Hence, 
and M are not affected by rotation of Y and Z; that is. 


'■a - V 
M = r 


( 31 a) 


Resolving equation ( 31 a) into components and using equations (ll) and 
(12) which define the aerodynamic coefficients give 


Cxi + Cyj + = Cy+i + Cy+j+ + 

c,T + Cml + c^k = cfl + cv+F + Cn+F 


( 31 b) 
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If j'*' and k"*" are resolved along Y and Z (eq. (30"b)) and the components 
are equated (eqs. (31I>))^ the following equations are obtained for the C' 
in terms of the C^*s and A. 


Ox 

0 

II 








Cy 

= Cy+ 

COB 

A 

- 

cz-^ 

sin 

A 



II 

0 

cos 

A 

+ 

c+ 

sin A 


Cl 

= Cz-*- 








Om 

II 

cos 

A 

- 

C ■’* 
'"n 

sin 

A 


On 

II 

0 

cos 

A 

+ 

Cn-^ 

sin 

A 

J 


(32) 


The components of the aerodynamic force_and moment along and 
are given by equations similar to equations (29) with all coefficients 
and all variables marked with the superscript +. Using this set of 
equations similar to equations (29) for the C+’s and substituting into 
this set the values of p'*', d'*'^ q+, and r+ given in terms of 

3^ a, 3^ d, q, r, and A by equations (30b), one obtains the following 
relations between the aenodynamic coefficients associated with the XYZ 
axes and with the XY^Z"^ axes from equations (32) 


= -CXo^ 






+ Zi 

q 2V 


-C 7 a - C 7 . — - C 7 
Ax 2V 2V 


cU ^ r^_ ^ 


Cl ^ 
^^p 2V 


-C7 ■‘a - C7»+ — - C7 -3d 
2V ^0. 2V 


= C 


+ 

Iv 2V 




(33) 


CniQa + Cm^ + Cmq 


p -t. p di . p + q^’ 

tni0, cx + OiH(^ -gy + Cniq 


pi 2^1 

-Cmcc,P - Cm,^ + Cmq — 


-Cm - Cm • ’*’ + Cm ZL 

oma P ^ma 2V ^ 2V 


J 


Since equations (33) must hold for all values of 3, ct, 3, a, p, q, 
and r, it is immediately evident that the aerodynamic coefficients asso- 
ciated with the XY'^’Z'*' axes, the C'*'’s, must be equal to the aerodynamic 



MCA TN 3350 


33 


coefficients associated -with, the XYZ axes^ the C’s. It may he con- 
cluded, therefore, that the aerodynamic coefficients are invariant with 
respect to rotation of coordinate axes about the axis of symmetry. Con- 
sequently, the Y and Z axes may be oriented at will around the X axis 
without regard to the orientation of the missile about X, since the aero- 
dynamic coefficients are solely functions of the missile’s external shape 
(and such nondimens lonal parameters as Reynolds and Mach numbers) and it 
has been shown that the aerodynamic coefficients do not change with orien- 
tation of the missile in roll. In fact, the missile may be allowed to 
spin about the X axis with respect to the YZ axes and the aerodynamic 
coefficients will be unaffected.^ 


Another consequence of the preceding analysis, which _is immediately 
evident, is that the aerodynamic force and moment, and M, are also 
invariant with respect to the orientation of the body in roll. The values 
of F^ and M change, of course, with spin through the rolling moment. 


C 


£i 

Zp 21 ’ 


and through the Magnus forces and moments, which will be discussed 


shortly. . 


A more adequate treatment of the consequences of rotational symmetry 
than has been presented here has been carried out in references 10, 11, 

12, and 13 . The general case of rotational symmetry about any submultiple 
of 360 ° is studied in these references in a most elegant and rigorous 
manner and it is shown that the full set of aerodynamic coefficients is 
reduced to the set given by equations (28) and ( 29 ) for rotational symmetry 
about all submultiples of 360 ° less than 180°, that is for an angle of 
360 '^/n with n = 3^ 5 • • • (any integer greater than 2). Ninety- 

degree rotational symmetry, or cruciform symmetry as it is commonly called, 
is only a special case and was chosen primarily to aid in visualizing the 
physical aspects of the analysis. Hence the phrase "or its equivalent" 
means rotational symmetry about any submultiple of 380 ° less than l80°. 

The 180 ° case lies in between the aircraft and missile cases but will not 
be treated in this paper (for an analysis of the l80° case, see refer- 
ence ih) . 


Geometry, Conditions of Flight, and Kinematic Relations 


The geometry is shown in figure 2. Three sets of Cartesian coordinate 
axes are used in developing the equations of motion, instead of two sets 
as formerly in the aircraft case. One set is fixed with respect to the 
earth and is known as before as "earth axes." The second set is fixed 


^The author is indebted to Mr. C. H. Murphy of the Ballistic Research 
Laboratories for pointing out to him the roll invariance of the aerody- 
namic coefficients of a roll symmetric projectile or missile. 
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partly In the missile and partly in space and. is known as "pseudo- 
Btahility axes." The third set is -fully fixed in the missile and is kno-wn 
as "body axes.” The relative orientation of _ .the three sets of axes is 
shown in figure 2 . Other features of the geometry such as the locations 
of the gravity and velocity vectors and the designations of quantities in 
the figures are the same as in the aircraft case and reference is made to 
the previous section "Geometry" and to figures l(b) and l(c). 

The earth axes are designated as before by XqYo^o- They are fixed 
with respect to the earthy but^ since the earth’s rotation is negligibly 
slow compared to the angular velocity of the missile, they are considered 
to be Galilean axes fixed in space. Their orientation in space is deter- 
mined by two factors: First, the Xq axis is alined parallel to the X 

axis in the neutral position of the missile; that is, the position for 
which the angles a, p, ijr, 0 , and cp are taken to be zero. In other 
words, the Xq axis is located by the initial conditions of the flight. 
Second, the Yq axis is horizontal and points to the right (of an 
observer facing forward) and Zq points down (but not necessarily verti- 
cally down) . As formerly, two locations are designated for the origin 
of coordinates of the earth axes , depending on the component of the motion 
being considered. For angular measurements, the origin of coordinates is 
taken to be the center of gravi-by of the aircraft at some particular point 
along its flight path. For linear measurements, the origin of coordinates 
is located as. required by the over-all scheme of measurement which would 
be set up to analyze the flight of the missile. 

The pseudo-stability axes are designated by XYZ. This set of axes 
exploits the freedom brought by the symmetry of the missile to orient the 
axes at will about the axis of rotational symmetry. The X axis lies 
along the axis of rotational symmetry. The Y axis lies in the space- 
fixed XqYq plane. Their origin of coordinates is at the center of 
gravity. In other words, the X axis is fixed in the missile and moves 
with it while the Y axis slides about in the XqYq plane. The positive 
directions along the axes are as follows; X is positive forward; Z is 
positive down; Y is positive in accordance with the right-hand screw 
rule. Should the X axis happen to be vertical in the neutral position 
of the missile, the orientations of the Y and Z axes (and the Yq and 
Zq axes) become arbitrary and would be determined by the conditions of 
the particular flight under consideration. It should be noted that these 
axes do not roll, since Y remains in the XqYq plane, and consequently 
the 9 and the derivatives of cp (dcf/dt, d^cp/dt^, etc.) of the XYZ 
axes are zero and remain zero d-uring the flight. It should also be noted 
that all three axes, X, Y, and Z, are principal inertia axes, as a con- 
sequence of the rotational symmetry of the missile. 

The body axes are designated by T, 2, 3. They are firmly fixed in 
the missile throughout the flight and are oriented so that 1 lies along 
the axis of rotational symmetry and ' 5 ' and "3 are coincident with Y and Z 
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at that point on the flight path which determines the initial conditions 
of the flight ( 1 ^ of course^ is coincident with X at all times). Their 
origin of coordinates is at the center of gravity. 

The orientation of the missile in space is defined by three angular 
coordinates^ two of which^ ^ and 6 , give the alinement of the pseudo- 
stability axes XYZ with respect to the earth axes Xq'^Zo; and one of 
whlch^ cp^ gives the alinement of the body axes 12 3 with respect to the 
pseudo-stability axes XYZ (see fig. 2). If these angles are all zero^ 
both the pseudo-stability axes and the body axes point in the direction 
of the earth axes. Any other orientation of the missile is reached by 
two consecutive rotations of the XYZ axes^ starting with XYZ coinci- 
^ent_with XqYoZq, and one rotation of the 123 axes^ starting with 
123 coincident with XYZ. In defining the rotations, the point of view 
is taken that the axes in guestion start in the position specified and 
proceed in ordered sequence to their final position. In particular, the 
two rotations of the XYZ axes are noncommutative and must be taken in 
.the order listed. 

ilf first rotation of XYZ. Start with XYZ pointing in the 

directions of XqYoZo. Rotate about Zq through t bringing 
Y to its final position. 

9 second rotation of XYZ. Rotate about Y through 0, bringing 
X and Z to their final positions. 

9 rotation of 12 3 « Start with 123 coincident with XYZ. 

Rotate about X through cp. 

The angles, 9 , cp, may also be defined as the angles between coor- 
dinate axes and the Intersections of certain planes . For the purpose of 
these definitions it is considered that the origin of coordinates of the 
XqYqZq axes is momentarily coincident with tha^ of_the XYZ axes . The 
point of view taken here is that the XYZ and 123 axes are in their 
final positions given by the orientation of the missile in space at the 
moment in guestion. The sign of the angle is specified by giving the 
axis about which the rotation is taken in going from the line named first 
to the line named second in the definition. The axis of rotation is 
listed in parentheses in each definition. 

angle between Xq and intersection of XqY^ and XZq planes; 
also angle between Yq and Y (rotation about Zq) 

9 angle between intersection of XqYq and XZq planes and X; 
also angle between Zq and Z (rotation about Y) 

cp angle between Y and 2; also angle between Z and 3 (rotation 
about X) 
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It should be. noted that the roll angle cp in the spinning-projectile 
and symmetrical -missile cases measures the rotation of the body axes 123 
and not the rotation of the pseudo-stability axes XYZ. As pointed out ' 

previously the XYZ axes do not roll. 

It may be of interest to note that the angles^ 6, are the 
Eulerian angles as defined in the classical treatment of the motion of a 
gyroscope (e.g., see section 43, "Heavy Symmetrical Top or Gyroscope," 
of ref. 15 ). The axis X corpesponds to the "axis of spin" (axis of 
rotational symmetry) of the gyroscope, and the axis Y corresponds to 
the "line of nodes." 

The following conditions are postulated for the missile and for its 
motion in flight: 

( 1 ) The missile has not only mirror symmetry, but also 90° rota- 
tional symmetry, or its equivalent, the axis of symmetry being the X 
axis . 

( 2 ) The magnitude of the velocity of ' the missile may be large but 
the change in the velocity must be small. Accordingly, the velocity may 
be represented by a constant plus a small perturbation; that is, 

V = Vq + u. (See eq. (2) .) 

( 3 ) The angles orienting the pseudo-stability axes with respect to ^ 

the earth axis, % 9, and locating the velocity vector, p, a, are all 

small. 

(4) The Y and Z components, q and r, of the angular velocity of 
the missile are small. 

( 5 ) The X component, p, of the angular velocity of the missile 
may be large but the change in p must be small. Accordingly, p may be 
represented by a constant plus a small perturbation; that is, 

P = Po + P* (34) 

(6) The flight path of the missile is a perturbation on a linear 
trajectory. 

Conditions ( 2 ), (3), (4), and ( 5 ) may be summarized by saying that 
the following variables must be small (e.g., have a numerical magnitude 
of 0.1 or less) and hence are first-order quantities; 



0 

if ' 



u 

P 

a 

^ small, for example, 
0.1 or less 

(35) 
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■where I is a characteristic length. Again it sho'uld be noted that if 
P, ql/SVo, and tI/ 2 Yq are first-order, then ^I/ 2 Vq and al/ 2 Vo will 
also be first-order. Also, it should be noted that Vq, cp, Pq, and 7 q 
may all be large. 

Again it should be emphasized that, although the linear theory gives 
an adequate description of the motion in flight of many projectiles and 
missiles, there may be cases in which the second-order terms neglected in 
the derivation of this theory are as large or larger than the first-order 
terms retained. Consequently, the analysis of an actual flight should 
include a check of the relative magnitudes of first- and second-order 
terms. The point here is that in both the aircraft and missile cases the 
a priori assumptions of the theory should be checked a posteriori by an 
analysis of the measurements . 

As in the aircraft case, the equations of motion are derived from the 
basic vector equations 



( 3 ) 

dt 

^ = M 
dt 

( 4 ) 


The components of the linear momentum, M, and the velocity, V, 
along' the XYZ axes are given as formerly by the tabular listing of 
equations (5)* 

In deriving relations for the components of the angular momentum, 

H, there are three angular velocities involved: 

(1) to angular velocity of the missile with respect to the earth 

axes, XqYqZq 

(2) ii angular velocity of the XYZ axes -with respect to the 

XoYqZo axes 

(3) w angular velocity of the missile with respect to the XYZ 

axes 

It can be shown that 


0) = IT + w (36) 

The components of are given by equations ( 7 b) with cp = (p = 0 , 
since the XYZ axes do not roll. For convenience in computing the 
relative magnitudes of terms in the equations of motion, these components. 
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may be written as follows, expanding sines and cosines and retaining 
only first- and second-order terms: 


Hy = r'QV 


= 


dt 




(37) 


Since the 1 and X 
given by 



j 

axes are coincident, the components of w are 



Wy = 0 > 

W2 = 0 


(38) 


Substituting equations (37) and (38) in equation (36), with 
gives the components of H and u along XYZ in terms of Q > 
by (correct to second order) 



I 

w 

i 

IXP 

p = ^ - 0r 

j 

lyq 

d0 
“ dt 

k 


r 

dt 


JX2 ” 0, 
and i|f 


(39) 


In equations (39) the transverse components of H, Hj and are written 
as lyq and I^r. How, it can be shown that rotational symmetry requires 
that ly = Consequently, a single value may be assigned for the 

transverse moment of inertia. This fact is used in deriving and solving 
the equations of motion for symmetrical missiles_j_ In deriving equations 
for p, a, 9, 0, and I'’ the rates of change of M and H will be computed 
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with respect to the XYZ axes and, hence, equations ( 3 ) and ( 4 ) must be 
transformed in order to account for the movement of the XYZ axes . The 
transformation for the rate of change of any vector A from fixed to 
moving axes is given in the missile case by 

= A + fl X A (40a) 

dt 

9 __ 

where A is the rate of change of A measured with respect to XYZ axes 
and is the angular velocity of the XYZ axes with respect to fixed 
(earth) axes. The components of the transformation are given by 


dt 


= Ax + qAg - rAy 


dAy • 

— Ay + rAx 


(40b) 


dt 


= Ag - qAx 


Aerodynamic ^d Gravitational Force Systems 


As in the aircraft case, the resultant external force F of equa- 
tion (4) is the sum of the aerodynamic and gravitational forces, as given 
by equation (lO), while uhe external moment M of equation ( 5 ) is solely 
an aerodynamic moment. 

Concerning the aerodynamic forces and moments, it is assximed, as_ 
formerly, (a) that the components of the resultant force and moment, 
and M, are given by equations (U) and (12); (b) that the coefficients, 
the C’s in equations (ll) and (12), are fimctions of the variables p, 
a, P, d, p, q, and r; and (c) that these functions may be expanded in a 
series in the variable named, where the general formula for any coeffi- 
cient, Ca (a = X, Y, Z, Z, m, n), is given by equation (13). Again, the 
coefficients of the series, the aerodynamic coefficients, are assumed to 
be independent of the dependent variables named and to be functions only 
of the body’s external contour and of such fundamental aerodynamic para- 
meters as Reynolds and Mach numbers. 

o 

In the aircraft case, all of the variables, p, a, p, d, p, q, and r, 
are assumed to be small, and consequently, all quadratic terms in the 
series expansion of the aerodynamic coefficients (eq.(l3)) second-order 

and may be neglected. However, in the missile case one of these variables, 
p, may be large and the quadratic terms of the series which Involve p 
may be first-order and must be included in the equations for the aero- 
dynamic forces and moments, equations (28) and ( 29 ). There are also 
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higher-order terms in p which are first-order in the other variables 
and therefore should be included from a strictly logical standpoint. How- 
ever, the available experimental data indicate that the forces and moments 
in question vary linearly with p to within the accuracy justified by a 
first-order theory. Accordingly, only the quadratic terms Involving p 
will be included. The general 
coefficient Ca^ which involve 
for X, Y, Z, Z, m, or n; 

Ca = ^ + 

dZ pZ 
2V 2V 


rZ pZ 
'-arp 2V 2V 

The quadratic terms involving 
as the Magnus forces and moments. 

The conditions of symmetry require that many of the Magnus coeffi- 
cients vanish. Furthermore, rotational symtoetry establishes relations 
between coefficients associated with motions in the pitch plane (XZ) and 
in the yaw plane (XY) . The Magnus terms remaining are listed below: 


formula for the quadratic terms of any 
p, is given as follows, where a stands 


^ + Ca- ^ + 


pzpz SiPl. 

+ ^app 2 v 2V 2V 2V 


(i^l) 


p are known in the ballistic nomenclature 


Cx 

- %P J 

Cy 

n vl' n dZ pZ qZ pZ 

^Yap^" 2V ^'^dp 2V 2V ^^qp 2V 2V 

Cz 

r ft Pr, P^ 

2V “ 2ipp 2V 2V ■ 2rp 2V 2V 

Cz 

(There are no rolling -moment Magnus terms) 

Cm 

p R P^ 0 . p • P^ 4 - p £i 

‘-mppP 2V “PP 2V 2V ^rp 2V 2V 

Cn 

pZ _ dZ pZ „ qZ pZ 

CnajP- 2Y + Cn<ip gy C^<3P pv 2V 


The equalities between coefficients are given by 
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Clap = Czpp 


"■ 

\ 


Cy" 

C]2 • ~ 


> 

(43) 

Cyqp = - ^Zrp 

Cnqp = - 





The aerodynamic forces and moments specified up to this point 
(eqs. (28), (29), ik2), and ( 43 )) arise from the interaction of the air 
with the missile’s principal aerodynamic surfaces, which surfaces are 
assumed to conform to the requirements of mirror and rotational symmetry. 
However, in practice there may be small asymmetries in the missile’s 
contour or asymmetries due to control- surface deflections, and under 
certain circumstances it may be desirable to include aerodynamic forces 
and moments caused by the asymmetries. Forces and moments of this nature 
will be omitted from the main development of the equations of motion 
carried out in this section. However, they are described in Appendix B 
and their effect on the motion is discussed briefly. For a more complete 
discussion of the consequences of aerodynamic asymmetry see reference 4. 

In summarization of the aerodynamic contributions to the external 
force and moment, the conventional aerodynamic coefficients are given by 
equations (28) and ( 29 ) and the Magnus coefficients by equations (42) and 

(43). 


The remaining contribution to the external force, the force of 
gravity, F^, may be resolved at once along the XIZ axes from equa- 
tions ( 17 ), since cp = 0 for the XIZ axes: 


"" “s(- sin 7 q - e cos 7^) 

F(jy = mg(i|r sin 7 q) > (44) 


^GZ = “is(- 0 sin 7 q + cos 7 q) 

.> 


Derivation of the Differential Equations of Motion 


The differential equations of motion for the projectile and missile 
cases are derived in a similar manner to the aircraft case. The kinematic 
relations and the equations' for the external force and moment are 
substituted into the basic vector equations ( 3 ) and ( 4) , ' resulting in 
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six scalar eq^uations, three associated with components of the force and 
three with components of the moment. Specifically, the substitutions 
involved are these: 

Equations (3) and (4) are first transformed from fixed to moving 
axes by equations (40b). Equation (3) for the linear momentum is devel- 
oped by equations (4) for the kinematic momentum relations and by equa- 
tions (10), (11), (28), (29), (42), (43), and (44) for the external force. 
Equation (4) for the angular momentum is developed by equations (39) 
the kinematic momentum relations and by equations (l2), (28), ( 29 ), (42), 
and (43) for the external moment. 

Again, the equation associated with the axial drag force (F^) Is 
treated as a special case and derived by resolving the vector equation (3) 
along the Xq axis rather than the X axis in order to avoid the 
(qA^; - I’Ay) term. The equations associated with the side force (Fy) and 
the normal lift force (F^) contain a term in mV which is given precisely 
by the axial drag equation. However, in the side-force and lift equa- 
tions, this term is approximated by -Cy pV^S/2, since it is believed that 
the remaining terms in the drag equation°are negligibly small under most 
circumstances. Full details of this step are given in Appendix A. 

The next two steps in the derivation are the same as in the air- 
craft case. The equations are linearized by separating into first- and 
second-order (and higher in order) terms and retaining only the first- 
order terms. The details of this step are given in Appendix A. Nondimen- 
sional time and concise notation are Introduced next. 

The final step in the derivation differs from that taken in the 
aircraft case in one respect and is similar in another respect. The neu- 
tral attitude of the missile is not the steady- state glide of the air- 
craft but is one in which the longitudinal axis of the missile points in 
the direction of motion (a = p = O) . Since symmetry requires the 

coefficient to be zero, the main lift force due to the effective angle of 
attack of the aerodynamic surfaces is lost in the neutral attitude, and 
the lift, drag, speed, weight, and flight-path elevation will not be in 
balance with the missile in this position. Consequently, the gravity 
terms cannot be replaced with corresponding aerodynamic terms as in the 
case of aircraft. On the other hand, it is appropriate to replace the 
factor, zct, by the lift and drag factors, ^oD^ utilizing thereby the 

static characteristics of the missile measiared in a wind tunnel. Symmetry 
suggests that it would also be appropriate to replace the factor y^ by 
the corresponding side-force and drag factors, ypy, ^oD? which would be 
measured in a tunnel. The reasons for making these substitutions are the 
same as in the aircraft case and reference is made to the discussion 
leading up to equations (l9)* The corresponding equations for the missile 
are 
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^ = ^OD 

+ Xqj, » 

^pL ~ ^OD 


(U5a) 


Substituting equation (28) into equation (45a) gives the relation between 
the side-force and lift factors of a symmetrical missile from wind-tunnel 
tests ^ namely 


ypL = -^oL f^ 5 h) 

The six differential equations are listed below in the same order 
as in the aircraft case: the three force equations axial drag, side 

force, and normal lift force, and the three moment equations, roll, pitch, 
and yaw. 


^ ^ ^ ^vxppDcp - v^xpp + xq = 0 (46) 

(1 - yp)Dp - ypj^p + (1 - yr)Dt - x^t - 

vydpDa - vycopa - vyqpD© =0 (4 t) 

(1 + Z(i)Da. + z^a - (1 - Zq)D0 + xq6 + 

vzppDP + vzppP + vzrpD^ - Zq = 0 (48) 

D^cp - ZpDcp - vlp = 0 (49) 

D^0 - mqD9 - m^J^Da - - vmppDp - 

vmppP + v(K - mpp)Dt = 0 (50) 

- Ui-Dt - npDp - npp - vn<ipDa - 

vnapCt - v(K + nqp)D0 = 0 (5l) 
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Discussion of the Differential Equations of Motion 


If equations (46) through (51) for the projectile and missile cases 
are compared with the corresponding equations (20) through (25) for the 
aircraft case, it is evident that they can no longer he divided into two 
separate and distinct groups in the same way that was possible in the 
aircraft case. Equations that formerly involved only the set of variables 
u/Vq, O', and 6 or the set cp, P, and ^ now contain members from both 
sets. Closer inspection discloses that the new members in the aircraft 
equations are all multiplied by the mean spin, v, and hence the spin is 
shown to be responsible for the interaction between that phase of the 
motion associated with w/Vq, o,, and Q and the phase associated with 9, 
P, and 'I''. In a word, the spin couples the longitudinal and lateral 
motions . 

Although the six equations can no longer be separated into two parts, 
one defining the longitudinal stability and the other the lateral sta- 
bility, it is possible to divide them into two distinct groups. Equa- 
tions (it6) and ( 49 ) involve only the variables u/Vq and 9 and constitute 
one group. Equations (4-7), (48), {^ 0 ) , and (51) Involve only the vari- 
ables p, a, 9, and 'k and constitute the other group. Furthermore, it 
will be shown that the later group can be reduced from foinr equations to 
two equations by introducing complex variables. Consequently, the differ- 
ential equations of motion for projectiles and missiles can be reduced to 
two distinct pairs of equations with each pair involving two dependent 
variables. The differential equations will be reformulated in this manner 
and will be solved in the following section. 

It may be instructive to rettirn at this point to a statement made in 
the introduction that "the aircraft equations do not describe the gyro- 
scopic nutation and precession of a spinning shell, and the projectile 
equations do not describe the phugold oscillation of aircraft" and to 
discuss briefly the reasons for this statement. The phugold oscillation 
of aircraft will be considered first. Equations (20), (22), and (24) may 
be simplified to give the essential features of the phugold as follows 
(see eqs. (lO-lOO) of ref. 5 ): 

D (1) » 0 (52) 

2-0 (f) - “ = ° (53) 

Comparison of equation (52) with equation (46) and of equation (53) with 
equation (48) shows that the z^6 and 2 zq(u/Vq) terms are missing from 
equations (46) and (48) for the missile and projectile case. Now, the 
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Zq factor is the constant lift coefficient that arises from the balance 
of weight by lift required by the equilibrium condition, of a steady- state 
glide. Since the equilibrium condition of projectiles and missiles is 
with their lift zero (or what corresponds to an equilibriim condition, 
since they are newer in a truly steady state in flight ifith all dependent 
variables zero, except in a terminal velocity vertical dive) , the Zq 
terms vanish from the projectile and missile equations and they do not 
predict the phugoid oscillation. In other words, projectiles are hurled 
rather than flown through the air. 

The gyroscopic nutation and precession of spin- stabilized projectiles 
will be considered next. The essentials of the gyroscopic motion are 
shown if the only external force or moment acting is taken to be the 
moment given by m^^, referred to in the ballistic nomenclature as the 
"overturning moment." In this case, the equations of motion reduce to 

B^6 + vKD'il'’ - ma/x - 0 (5^) 

, vKDe - npp = 0 (,55) 

Comparing equation (.5^) 'with equation (24) and equation (55) with 
equation ( 25 ), shows that the vKDi|'' and vKDe terms are missing from the 
aircraft equations. Since these particular terms are essential to the 
gyroscopic motion (as will be shown shortly), it is clear that the air- 
craft equations cannot predict the motions of rapidly spinning projectiles. 

Dynamic stability requirements for spin-stabilized projectiles are 
not as widely known among aeronautical engineers as those for aircraft 
and it may be informative to discuss this aspect of the subject briefly. 
The essential features of the gyroscopic precession and nutation can be 
developed from equations ( 5 ^) (55) • this simplified case, the 

traJ_ectory is a straight line, since there are no external forces^ X, Xq, 
and V are all coincident j and the angles of pitch and yaw are related 
to the angles of attack and sideslip as follows: 

P = -t 
a = e 

Using equations ( 56 ) and the equality np = -mo, from equations (28), 
equations ( 5 ^) (55) become 


D^a - vEDp 

- iiiqP' = 0 

(57) 

D^P + vKDa 

- niccP = 0 

( 58 ) 
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Equations (57) aii-d (58) can be reduced to a single equation in the 
complex variable, 

I = P + ia 

namely, 

- ivKDl - ma,^ = 0 (59) 

where i(WrT) has the significance of an accounting parameter permitting 
the simultaneous solution of the pitch and yaw equations . 

The solution of equation (59) is given by 


g = A^e + Age (60) 

where Aj_ and A2 ai’e constants 
ditlons and of and S^} and 

ivK + V- + 4m„ 

S = a 

^ 2 

It should be noted that met, is positive for spin-stabilized projectiles 
(i.e., the center of pressure is ahead of the center of gravity; if met 
were negative, the projectile would be arro^v stable and spin would not 
be needed for stability) and the radical will be either real or imaginary 
depending on whether is smaller or greater than 4ma,« The other 

factor in equation (61), ivK, is always imaginary and hence the motion is 
always an oscillation. 

Now, if the radical in equation (61) is real, the S’s will be 
complex and one of them will have a positive real part. Hence, the 
oscillatory motion will diverge. On the other hand, if the radical is 
imaginary, the S*s will be purely imaginary, and the motion will be 
an oscillation at constant amplitude, which is the stable type of motion 
in this simplified case since damping terms have not been included. Con- 
sequently, the criterion for stability is 


which are functions of the initial con- 


St and 


S, = 


are given by 
ivK -nZ-K^v^ + 


Hi 


(61) 


>^a 


or 


4ma, 


> 1 


In ballistic nomenclature, is called the "stability factor." 

In other words, the analysis of the simplified equations of motion has 
led to the well-known requirement for stable spinning shell, namely, that 
the stability factor must be greater than unity. 
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The complete differential equations of motion are solved in the 
following section of this report. In the solution, the equations are 
transformed from time to distance as the independent variable. However, 
it can be shown that the nondimens ional equations in t"^ are the same 
as those in x+ except for minor differences in coefficients, which will 
be accounted for in the development that follows . 

The integrated equation for the I (f^) history will have precisely 
the same form as equation .( 93 ) for the SCx"*") history. It can be seen 
that the complete | equation has the same form as the simplified equa- 
tion (60) except for an additional term, the constant, Hq. For the 
^(x"^) equation ( 93 ) "the relations between the exponents in the equation 
and the aerodynamic factors are given by equations (l 62 ) and (163). For 
the |(t+) equation (60) the corresponding relations are given by 


®iR “ + ^ 2 .} “ 'I'l “ 


'll 


vK(l+a) o vK(l - 0-) 

- = . - ! 


(62) 


where 


^1 = 




= 


2 m 


2 q 


PP 


-/ 


1 - 


s = 


Wa 


= stability factor 


The above solution 
is valid for -co<s<0, 
that is, for missiles 
with arrow stability cor- 
responding to ma being 
negative, and for 
Ks<+co, that is, for 
spin-stabilized projec- 
tiles with a stability 
greater than unity. The 
values of cr correspond- 
ing to values of s are 
shown in sketch (c). The 
solution is not valid. 



0.0 1.0 

Sketch (c) 
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however^ in the region 0 < s < 1 (shown crosshatched in the sketch) for 
which a is imaginary. 

Equations (62) are in the most suitable form for spin-stabilized 
projectiles (s > l) . However, they have a clearer physical significance 
for missiles with arrow stability (mo, < O) if they are rewritten as 
follows : 


Sir = + T2 


SaR - - T- 


Sil = I (vK + ^-kmc,+ v2k2 ) , 


(63) 


^2l 


I ( vK - 


where 


Ti = 


(^L " “q “ 
2 




vK 

2 



+ mq 




K 


) 


Equations (62) show that the criteria for the dynamic stability of 
spin-stabilized projectiles, that is, the requirements for a convergent, 
oscillation, are 


> 1 (or a > 0) 




SgR < 0 


(6it) 


The term Tx is normally negative. Hence, it can be seen that dynamic 
stability requires that T 2 < Tx regardless of the sign of Tg. lu 
practice, the magnitude of Tg depends on a balance between the lift, 
z^j j the damping moment, mq + mo,; and the Magnus moment, m^p. That spin- 
stabilized artillery projectiles are in reality fully dynamically stable 
is a fact that has been thorou^ly verified by experiment. 

Lift and Magnus moment are both vital elements and must be included 
in the equations of motion in order to derive the correct criteria for 
dynamic stability. If they are omitted, equations (62) take the form 
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Now 0 < cr < 1 for 


(mq + m^) 
" 2 

(mq + m^ 
S2R = -^2 

1 < s < 00 ; hence. 





- 

, 


(65) 


Therefore, one S will be positive even though (mq + m^) is negative, 
and a simplified analysis with lift and Magnus moment omitted would lead 
to the erroneous conclusion that spin-stabilized artillery projectiles 
cannot have a convergent oscillatory motion (see section on "Application 
of Results: Missiles" of reference l4) . 


APPLICATION TO THE FREE-FLIGHT TEST FACILITY: 
THE AERODYNAMIC RANGE 


During the past two decades a new facility, the aerodynamics range, 
has been developed for measuring the aerodynamic properties of bodies in 
free flight. The range is properly classed as a flight-test facility, 
since the aerodynamic measurements are made during a completely free 
flight of the model. Accurate records are taken of the model’s movements 
along a certain length of its flight path and the aerodynamic characteris- 
tics are determined from these records . 

The experiment consists of recording the positions, angular orienta- 
tions, and times of the model at a series of stations placed along its 
flight path through the range . Photography is the primary medium used for 
recording, since it is a precise technique and one which does not inter- 
fere in any way with the model’s flight. An electrical spark discharge 
generates the light for the photography and its duration can be made so 
short that the picture is nearly instantaneous in relation to the movement 
of- the model during the time of exposure. Hence the photographic record 
gives the x, y, z, cp, 0, and of the model at the particular Instant 
of time that the spark produces the exposure. The remaining flight datum, 
t, the time at which the photograph is taken, is measured by a special 
high-precision chronograph. To enumerate specifically, the experimental 
data from the flight records are y, z, <p, 0, and f at a series of x 
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and t along the trajectory. In addition^ measurements are made of the 
model’s physical properties^ that Is, its dimensions, weight, center of 
gravity, and moments of inertia, before flying it in the range. 

To the author’s knowledge, there are seven aerodynamics ranges in 
operation at the present time: The Aerodynamics and Transonic Ranges and 

the Controlled Temperature-Pressure Chamber at the Ballistic Research 
Laboratories (U. S. Army), the Pressurized and Aerodynamics Ranges at the 
U. S. Naval Ordnance Laboratory (U. S. Navy), the Aeroballlstlcs Field 
Laboratory at the U. S. Naval Ordnance Test Station (U. S. Navy), and the 
Supersonic Free-Flight Wind Tunnel^ at the Ames Aeronautical Laboratory 
(MCA) . The forerunner of all these facilities, the Aerodynamics Range 
at the National Physical Laboratory, Teddington, England, is dismantled 
at the moment of this writing, although it is understood that plane are 
in effect for its reconstruction at Fort Halstead (Kent, England). 

The aerodynamics range is unique among fll^t-test facilities not 
only in its measurement techniques but also in the conditions under which 
its testing is carried out. The extent of the flight path is severely 
limited . The region of space tinder observation varies from a length of 
15 feet with a cross section 1 -foot square in the smallest range to a 
length of 750 feet with a cross section 25-feet square in the largest 
range. Testing is confined to flights for which the trajectory is nearly 
a straight line and the changes in velocity and angular inclination of 
the model over the length of the range are small. Most of the models 
themselves are either simple bodies of revolution or bodies with cruciform 
wings and fins having 90 ° rotational symmetry. Consequently, test condi- 
tions in the range agree with the flight conditions assumed for projectiles 
and missiles in the present analysis, and the equations derived herein 
describe correctly the linear and angular motions of models flown through 
the range . 

Two steps are required to obtain the aerodynamic characteristics of 
the model from the flight-test data. First, the constants in the equa- 
tions of motion are evaluated to give the best possible "fit" to the 
experimental measurements. Second, the aerodynamic coefficients are com- 
puted from certain of these constants. In this section of the report, the 
differential equations of motion will be solved for the particular test 
conditions prevailing in the aerodynamics range, the process of "fitting” 
the equations to the test data will be discussed briefly, and relations 
between the constants of the equations and the aerodynamic coefficients 
of the model will be derived. 


2 

At the Ames Laboratory, the test chamber of the range is the work- 
ing section of a supersonic wind tunnel. Thus the Ames facility combines 
the wind stream of a supersonic wind tunnel with the measurement techniques 
of the range. The purpose is to carry out free-fllght testing at much 
higher Mach numbers than would be possible in still air at normal room 
temperatures (see ref. I6). 
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Equations of Motion 


The conditions of flight are modified slightly to fit the particular 
circumstances of testing in the range. The trajectory is so nearly hori- 
zontal that 7 o wilL be considered equal to zero in the Xq and zq fac- 
tors^ thereby making 


X(j = 0 




mg 

pSVo^ 


( 66 ) 


The variation of drag with spin is so small according to the limited 
experimental evidence available (see ref. I7) that the xpp terms may 
be neglected in the drag-force equation, (^ 6 ), compared to the Xq term. 
However, should the conditions of any particular test be such that the 
Xq and Xpp terms are not truly negligible, they may be readily included 
in the general solution of the equations of motion since both are linear. 


The independent variable in the equations of motion will be changed 
from time to distance and the equations transformed accordingly. It has 
been claimed that distance is the natural variable rather than time (see 
ref. 11). From the experimental standpoint, the distance is recorded in 
every photograph while time is recorded only at a few, widely spaced sta- 
tions along the range . From the theoretical standpoint , the Inertial 
forces vary with the square of the time and, in transformation to distance 
as independent variable, with the square of the velocity. Now, the aero- 
dynamic forces also vary with the squane of the velocity, so it might be 
expected that the velocity would disappear as a primary parameter from 
the transformed equations . 

These advantages are realized in the case of the rolling motion 
(eq. (ij- 9 ) for 9), for which distance is the better variable on which to 
base the analysis than time. The velocity no longer appears explicitly 
in the <P(x) equation and the speed range over which the equation correctly 
describes the rolling motion is limited only by the variation of the aero- 
dynamic coefficients with Mach and Reynolds numbers. Furthermore, the 
basic differential equation for the rolling moment is linear in cp and 
contains no first-order term involving the other dependent variables . It 
is not necessary to limit the variation in the value of 9 in order to 
linearize the equation and its solution is valid for large changes in 9 
as well as for small changes . Also the rolling motion is Independent of 
the other motions and may be treated as a special case. Accordingly, a 
constant aerodynamic rolling moment, is introduced in order to 

include the class of missiles having aerodynamic surfaces with aileron 
deflection. Of course, aileron deflection violates the requirement of 
mirror symmetry but is permissible in the rolling-moment equation since 
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relaxing this requirement introduces no nev terms in the rolling-moment 
equation. It is to be understood in the subsequent analysts that only 
the cp(x) equation is valid for large changes in spin and possibly also 
in velocity; all other equations are still subject to the limitation of 
small changes in both spin and velocity. 

Actually, the advantages of using distance rather than time are not 
as great as might be supposed. The experimental data are reduced in such 
a manner that time can be computed easily and accurately for all of the 
photographic records. The transformed equations, vith the exception of 
the rolling-moment equation (^4-9) for cp , still contain terms Involving 
the velocity, and it is necessary to postulate that the change in velocity 
be small in order to linearize the equations of motion. In fact, the 
equations for the transverse displacement and for the pitching and yawing 
motions are practically Identical for either the nondimens ional time or 
the nondimens ional distance as independent variable, and flight data from 
the range could be analysed on either basis. However, in conformity with 
accepted ballistic practice, distance will be selected as the variable 
used in developing the equations of motion for the range. 

In developing the transformation from time to distance it will be 
recalled that the x distance is measured along the space-fixed axis, 

Xq* On the other hand, the velocity is strictly determined by the rate 
of change of distance along the trajectory, that is, letting s be the 
arc length along the trajectory 

V = ^ 
dt 

Hehce, the relation between V and x is given by 



where 

If = sec < VXo = 1 


+ .-V + • • • 


Now, the angle between V and Xq is normally very small, being less than 
2° ( 0.03 radians). Therefore the term (^VXq)^/2 is of the order of 
magnitude of 0.001 or less and may be safely neglected compared to unity. 
Hence, the velocity may be approximated by 


V 


dx 

dt 


( 67 ) 


Using equation (67) for the velocity gives the transformation from time 
to distance for any dependent variable A 
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dA _ y M 
dt dx 

d^A yg d^A dV dA 

dt^ dt dx 


( 68 ) 


Wow dV/dt is given by equation (A 18 )^ (since the Xpp and Xn tenns 
are neglected) as follows 


dt 


m 




( 69 ) 


After equation (69) is substituted in equation (68) the transformation 
becomes 


M ^ y M 

dt dx 

► 

^ = V2 ^ yZ Ml 

dt^ dx^ “1 


(70) 


Distances are nondimens ionali zed in a manner similar to times by 
defining a base distance, \iZ , and denoting the nondimens ional distance 
by a plus; that is, 


v+ = 


X 

p.1 


7 ^ = 


X 




(71) 


Differentiation with respect to 
quantity, A 


z"*' = 


pZ J 

is denoted by 


ID> ; that is , for any 


EDA 


dA 

dx"*" 


(72) 


Hence, the nondimensional form of the transformation from time to distance 
is given by 


DA = 



©A 


D^A = 



( 73 ) 
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One possible procedure at this point would be to substitute equa- 
tion (73) into equations (46) through . (51) ^ thereby transforming directly 
from t'*' to x+. However; it can be seen that neither (v/Vq) nor (v/Vq)^ 
is a common factor in any of the equations of motion; (46) through (51) ^ 
and that; if this procedure is followed; it will be necessary to make 
further approximations in order to linearize the equations. A better 
procedure is to return to the exact vector equations; (3) and (4); the 
kinematic relations; equations ( 39 )> and the moving-axes transformation; 
equations (40); to transform these from t to x using equations (?0); 
and to derive the noridimensional linearized equations in x"*" as formerly 
in f*". The detailed steps are carried out in Appendix A. 

Two new equations in addition to the set corresponding to equa- 
tions (46) through (51) are needed for the analysis of flight tests in the 
aerodynamics range. The trajectory is recorded and serves to determine 
the aerodynamic lift coefficient of the model. The position of the model 
transverse to the Xq axis is given by the y and z coordinates of its 
center of gravity measured along the Yq and Zq axes. In the derivation 
the force system is simplified by retaining only the lift; drag; and 
Magnus aerodynamic forces together with the gravity force; since experi- 
ence in reducing experimental data from the range has shown that the 
contributions to the transverse motions by the remaining components of the 
aerodynamic force are less than the errors of measurement and; hence; may 
be neglected. The detailed derivation of the y and z equations is given 
in Appendix A. 

The eight differential equations of motion are listed below with the 
first three being the equations for the trajectory (x; y, and z) and the 
remaining five for the transverse velocity (3 and a) and for the angular 
motion ( cp; e ; and i(r) , It will be recalled that rotational symmetry estab- 
lished certain equalities between aerodynamic coefficients associated with 
motions in the pitch and yaw planes and; taking advantage of this fact; 
the notation will be simplified by using only the aerodynamic coefficients 
associated with the pitch plane. The selection is arbitrary; of course; 
and the yaw plane coefficients could have been chosen Just as well. 


D^X+ + Xq (1 + 

Xot+)"^ = 0 

(74) 

+ Zqj^p - 

vzppu » 0 

(75) 

P^z+ + + 

VZppP - Zq 

(76) 


(1 + Zct)lD)p + ZqlP + (1 - zq)lD)t - vz^Pa - vz^^a + vz^pPe = 0 


( 77 ) 
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(1 + Z(5jPa + Z^a - (1 - Z(j)E2)e + vz^pDp + + vz^^Dijr = Zq (78) 

nj^cp - (Zp + Xotocp = lo (79) 

ED^0 - (mq + Xq) EDe - ct - - vm^pOP - ^TUppP + v(K - mpp)Pi(f = 0 

( 80 ) 

2 

© t - (mq + Xo)©t + m^P + ma_p - vmppIIDa - vm^pa - v(K - ©e = 0 

( 81 ) 


Inspection of equations ( 75 ) and (76) shows that the coefficients^ 
^L-* ^^PP ©^y^^ Pj and a, respectively, in equation ( 75 ) 

reappear as the coefficients of ID)^z+, a, and p in equation (76). The 
same similarity of coefficients occurs between equations ( 77 ) and (78) 
and between equations (80) and (81) and suggests that the three pairs of 
equations may be reduced to three single equations by a proper choice of 
dependent variables. This, in fact, is the case; by defining three new 
dependent variables 


A = y"*" + iz+ 


I = P + la 


) (82) 


T] = ilr - le (ip = 0 + i\)f) j. 

and including equations (7^) and (79) in order to complete the set of 
differential equations of motion, equations (7^) through (8I) reduce to 


T -.2 4 - 

D x^ 


+ Xr 


(1 


+ X^f^) 


= 0 


(83) 


ID)^A + z ^5 = izg 


(8h) 


z^ID)fl + z.DJI + z^^ = Izq 

(85) 

ED^cp- (ip + Xq )©9 = Iq 

( 86 ) 

©^p - m^©P + m»©^ ■+ m^l =0 

(87) 
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where 



“ ^aL 


z • 

= 1 + 

z. -f ivz, 

1 


a 1 

Zti 

= 1 - 

Zq + ivz- 

4JL/I 

= m_ + 

CL 


m. 

= m. + 

ivm. 

1 

a. 

^P 

mr| 

= mq_ + 

Xq + iV 


rp 





( 88 ) 


The physical significance of the complex dependent variables^ ^ ^ 
and T]^ is evident if they are considered to be vector quantities. The 
quantity A is the vector displacement of the trajectory transverse to 
ttie Xq axis . The quantity | is the component of the -vector velocity^ 
Vj transverse to the model's axis, X; in ballistic nomenclature | is 
known as the "vector yaw." The quantity Hi is the component of a unit 
vector along the model's X axis transverse to the space-fixed Xq axiSj 
thereby defining the orientation of X axis in space. Because of its 
rotational symmetry^ the model's motion in flight can be described by 
these three vector quantities together with the roll angle, cp; 
time-distance history, f^(x"*’). 


The solution of the differential equations of motion is quite 
straightf onvard . The coefficients in all equations are constants. Equa- 
tion (83) contains only the dependent variable x"*" and, although not 
linear, may be Integrated in closed form. Equation (8^4-) is linear in the 
dependent variables A and i and may be integrated once a solution has 
been obtained for i. Equation (86) contains only the dependent variable, 
cp, is linear in cp, and may be integrated directly. Equations (85) and 
(87) form a pair of simultaneous, linear differential equations for | 
and T) and may be solved by a variety of methods; for example, by the 
use of the Laplace transform (see ref. I8). 

The exact solution to the x"''(t'*') equation (83) is given by 

x+ = In (l+Xot'*') (89) 

The logarithm in equation (89) is expanded in a series, since is 

ordinarily a small quantity, and x is expressed as a function of t 
rather than x+ as a function of t'*' so that Vq can be determined 
directly from the time-distance records as is required for the computation 
of T. Since x+ = 0 , V = Vq, cp = cp^, ©cp = v, A = A^, ©A = (©A)^, 

I = Ti = t]q, and ©Tj = (©ti)q at f*" = 0, the solutions to the differ- 

ential equations of motion (eqs. 83 through 87) are given as follows; 
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X = Sq + + a^t + agt^ 


A = + b^x"^ + bgCx'*')^ + bg 

, CgX 

9 = Gq + 

. , 3 dox"^ , 

I = do + dj_e + dge ^ 


d^2 * dt ,^ 


doX"* 


T] = fg + f^x"^ + fa® ^ + ^*3® 


,(34 x"* 


5T 

(90) 

(91) 

(92) 

(93) 

(9^) 


where 

ag = constant required to adjust the zero of the distance scale to coin- 
cide with the zero of the time scale ^ 


X = 0 at t = 0 
a = V = 

- Vg - ^ 

a, = - 

aa 2 t 

X sy 

a = ■ ■ - Q- -- Q ' (assuming Xg constant) 
3 t -2 


_ / It ( 3g 

tlQ = 

dg d4 


b^. = (PA)g 4 Zg^|i+ ^ 


bs 

bo 


izQ - Z|dp 


®o = 9o - 


v(Zp + Xg) + I-c 

(ip + Xq)^ 


-Zf 


(Zp + Xg) 


(95) 

(96) 

(97) 

(98) 

(99) 

(100) 

(101) 

(102) 

(103) 


(loh) 
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Co — 


Co = 


v( Zp + Xq) + Z( 
(Zp + Xq)^ 




dn = - 


d. = 


d.^ = 


izgHiTi 

dgd 4 

P(d 2 ) 


*^2(^2 “ ^^4) 
P(d 4 ) 

d4(d.4 - d. 2 ) 


P(S) 


+ 


iz. 


S = 


Zi " 

I ^ 

d 2 ;d 4 

roots of_ R(s) = 0 


:t) [(IG)ii)o + ] 


z » 

i 


VS - 




*^2 ^*^4 

r(s) = (S - dgjlS - d4) = 0 
r(s) 


= 32 + 


i’o = - 


Z " HIyi 2 • “ 

I n I T1 I 


Q(dp) 


S - 

Q(d. ) 


m^z^ + 

z 7 


d^^ (d2 - d^) d/ (d^ - d„) 


iz, 


f, = 


fo = 


G®| 


z^dad 4 

Q(da) 


2 4 2 


V (^2 - ^4) 


Q(d 4 ) 
‘3 d^2 (3^ . 


f = 


(105) 

(106) 
(lOT 

(108) 

(109) 

(no) 

(in) 

(112) 

(113) 

(114) 

(115) 

(116) 

( 11 ?) 


(118) 



KACA ra 3350 


59 


Q(s) = (Ptj)o S2 + 


(©Tj)oZt + m,2 5 - m 


g^gSp - mgZ|io 




lz„in 


"I 


S - 




(119) 


Eguations (90) through ( 9 ^) with the constants defined by equa- 
tions (95) through (119) describe completely the motion of a projectile 
or missile flying through the range under the conditions postulated. They 
form the basis for determining the aerodynamic properties of the model 
from an analysis of the flight records . 


Reduction of Flight Data 


The first step in the reduction of data is the computation of the 
constants appearing in the equations of motion. It is assumed that the 
equations represent correctly the actual motion of the model and that any 
differences between the theoretical and experimental values are due to 
errors of measurement. The problem is to determine the values of the con- 
stants, the a'sj b’s, c's, d’s, and f’s, that give the best fit of theory 
and experiment. Specifically, the best fit normally means that the sum 
of the squares of the residuals between the computed and measured values 
is a minimum. 

In principle the procedure is to compute an Initial set of values of 
the constants by some approximate analysis of the data and to correct this 
initial set through a series of iterations. The sum of the squares of 
the residuals is computed at each step of the iteration, and the varia- 
tion of the residuals with time or distance is studied in order to detect 
any systematic trends. If the process is convergent, it is carried out 
until the residuals are random in their distribution and until the sum of 
their squares approaches a steady, minimum value. Then, if the average 
value of the residuals after the final iteration is of the same order of 
magnitude as the estimated experimental error, the fit is considered to 
be satisfactory. 

It should not be supposed that any one set of procedures has been 
drawn up to reduce the data from aerodynamic ranges . The computations are 
detailed and lengthy, and each range has its own unique methods adapted 
to the peculiarities of the apparatus and to the uses for which the data 
are obtained . Reference is best made to the facility of Interest for a 
description of the methods in use there. 
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The final step in the reduction of data is the computation of the 
model's aerodynamic coefficients from the constants of its motions. The 
relations between the two are contained in equations ( 95 ) through (ll 9 ). 
However^ the equations named are formulated from the standpoint of com- 
puting the motions after the aerodynamic properties and the Initial con- 
ditions have been given. It is desirable to transform these equations 
and to simplify certain of them in order to facilitate the computations 
involved in this step. 

With regard to the constants in the x(t) equation ( 9 ^), the a*s^ 
the velocity and time factor are given by: 

Vo = (120) 


T 



and the drag coefficient factor by: 


Xo = 







( 121 ) 


( 122 ) 


With regard to the constants in the A(x"*') equation ( 91 )^ tiie b’s^ 
the lift coefficient factor^ and Magnus force coefficient factor are given 
by the real and imaginary parts of bg as follows: 


^oL " “ ^3R 


(123) 


= - 1 > 3 I 


( 124 ) 


although it should be noted that the determination of Z|jp from bgj 
may be marginal under certain circumstances. 

With regard to the constants in the 9 (x'^) equation ( 92 ); the c's, 
the static and damping rolling-moment coefficient factors are given by 

Ip + Xq = cg (eq. (106)) (125) 


^o “ 


and the initial value of the nondlmensional spin by 


V = Cj_ + c^c, 


z 3 


( 126 ) 


(127) 
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With regard to the constants in the equation ( 93 )^ "the d' b , 

the relations between these and the aerodynamic coefficient factors are 
derived by equating the sum of the roots of the quadratic in equa- 
tion (113)^ to minus the coefficient of the linear term^ equation (ll 4 )^ 
and the product of the roots to the constant term; that is , 


z.- (d, -^dj 

(128) 



^ 

Z; 2 4 

( 129 ) 


Substituting equations (88) into equations (l 28 ) and (l 29 )^ regroup- 
ing into real and imaginary parts ^ and replacing the aerodynamic factors 
by their definitions in terms of the aerodynamic coefficients gives for 
equations (l28) and (l29) 



- - (‘3aR + + <i4i) 


(130) 
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{ 


_2Ky2 





( 




+ 



1/2 


8\i^ 


4Ky^ 


V 


2 





+ iv 



■^m 


PP 


4Ky2 


SiiKy^ 2 


4^1 \ 4Ky2 2 J 


4|a 4 Ky 2 8 (i 


•'ffl, 


a 


2 Ky 2 



Cz 




<^4R 



+ i 



‘^4R + <^2R '^41 


( 131 ) 


It can be seen that many of the terms in equations ( 13 O) and ( 131 ) 
contain the projectile density factor^ [i. Now^ p. is a large number in 
all practical cases. It varies from 100 to 1^000 for full-scale projec- 
tiles and missiles and from 1^000 to 10^000 for models tested in the 
range. Consequently, certain terms will be much smaller than other terms 
and may be properly neglected consistent with the other approximations of 
the linearized equations. 

Separating equation ( 130 ) into real and imaginary component equations 
and neglecting all terms of magnitude l/(i or smaller gives 


(zqL - Xq - fflq - m^) = - (^2R + <^4 r) 


( 132 ) 


vK - (dgi + 


( 133 ) 


Separating equation ( 131 ) into real and imaginary component equa- 
tions , retaining only the terms of magnitude p in the real part , neg- 
lecting all terms of magnitude l/p or smaller in the imaginary part, 
and neglecting dgR d^R compared to dgp d^j gives 
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“a = <^4l (13^^) 

®Pp^ ^ (*^21 ‘^ 4 R + fisR ^541) (135) 

The sum of the damping-moment factors ^ mg + m^^ can be computed from 
equation (l32) using equations (l22) and (l23J for Xq and z^. The mean 
nondimensional spin v can be computed from equation (l 33 ) and checked 
against the value from equation (l27). The static-moment factor^ mo;,^ is 
given by equation (13^) • The Magnus moment factor ^ iJlpp> can be computed 
from equation (l 35 )* 


From equations (7^) through ( 81 ), it can be seen that all of the 
aerodynamic coefficient factors have been accounted for thus far in the 


reduction of data except 


z . , 

a’ 


'Q.^ possibly, 

indi- 


m^, can be computed, but not 


mo and m. 
y. a 


Also, only the sum, m^ + 

vidually. These exceptions represent the limits of the analysis of test 
data presented up to this point. 


One further possibility may be exploited in flight testing in the 
aerodynamics range. Instead of only one model of a given configuration 
being tested, two models are tested, both having identical external con- 
tours but one having a center-of-gravity position different from the 
other. In this way the sum of the damping-force factors, z. + Zq, may be 
determined from measurements of the sum of the corresponding damping- 
moment factors, + nig, at the two center-of-gravity positions; similarly, 
the Magnus force factor, may be determined from measurements of the 

Magnus moment factors, at the two center-of-gravity positions. Also, 

the normal-force factor, may be determined in this way from the ma’s, 
should measurements of the trajectory be unsuitable for some reason. 


On the other hand, it does not appear possible to separate + Zg) 
oi" ®q) or to determine ^rp-’ ™rp range tests. 

However, it can be seen from the development leading from equation (l28) 
through to equation (l 35 ) that only the sum (z. + Zq) or the shm (m. + mn) 
is retained in the final equations as a quantity having a significant 
magnitude and that all of the factors, 2 ^^, z^p, J“^p^ ®rp ai'® negli- 
gible. Consequently, it may be concluded that all of the aerodynamic 
coefficients of practical significance may be determined from free-fll^t 
testing in the aerodynamics range on projectiles and missiles with rota- 
tional and mirroi symmetry flying under the conditions of small angles of 
sideslip, attack, pitch, and yaw, and small changes in velocity and spin, 
as specified. 
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Comparison With Ballistic Theory 


Differences in nomenclature and derivation complicate a comparison 
of the projectile equations developed in this paper with those of ballis- 
tic theory. The differential equations giving the x(t), cp(x+)^ and ACx'^) 
histories can be compared directly with the corresponding ballistic 
differential equations^ as in the aircraft case. On the other hand, the 
equation giving the l(x'*') history is best compared in its integrated 
form. It is first necessary to derive the relations between the nomen- 
clature of this report and that used in ballistics . The comparison will 
be made by transforming the equations of this report into the ballistic 
terminology. The .ballistic nomenclature will be taken from reference 10, 
except for that used in the cpCx"*") equation, which will be taken from 
reference I 9 . The ballistic equations will be taken from the following 
sources: the x(t) equation from reference 20; the cp(x"^) equation from 
reference 19 ; the A(x+) and |(x+) equations from reference 21. These 
references were selected because it is understood that they constitute 
the basis for the analysis of data from the Aerodynamics and Transonic 
Ranges at the Ballistic Research Laboratories and hence present a section 
of ballistic theory which is in use at the present time. 


In ballistics® the aerodynamic coefficients are denoted by the capi- 
tal letter K and the aerodynamic factors by the capital letter J, the 
two being related in all cases by 



(136) 


where d is the diameter of the projectile. A study of the geometry and 
nomenclature of reference 10 shows that the following relations hold for 
the transformation of quantities defined in this report into ballistic 



®If a ballistic symbol is not specifically defined in this section, 
it has the same definition in both systems of nomenclature and is given 
in the list of symbols at the beginning of this report. 
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The transf oi'ma.'tion of the aei'Od3Tnaiiilc factors can be derived from the 
definitions of the ballistic factors^ the J’s, given in reference 10^ 
namely 


T7 - T 






?2 + iFa = ^ (- Jn + ivJp)(u2 + iug) + mu(vJj^ + iJg)(t02 + iug) 


G = -mdu 


>(138) 


Bu 


G2 + iGg = -' g - ■ (- vJqi - iJj4)(u2 + iug) + 
k d 


Jg + iVjYTn)(c02 + iCOo) 


k^d 




If equations (138) for Fj_^2^3 ^1^2 3 equations (ll), (l 2 ), ( 28 ), 

(29), and {k 2 ) for az %,Y,Z substituted into the rela- 

tions between the force and moment components given by equations (l 37 ) an <3 
the coefficients of the dependent variables are equated, since equa- 
tions (137) aiust hold for all values of the dependent variables, the fol- 
lowing relations are obtained: 


^oD " 

% = I" 


zq 


+ z. 
a 


Js 


^pp = - Jf 
Zrp = “ 11(1^ 



+ m. 
a 






> 


(139) 


^ ^ T 



It will be noted that the factor z. has been added to Zg in the equa- 
tion for Jg and similarly the factor m^ has been added to mg in the 
equation for Jq. In the strictest sense the added terms should have 
been omitted since the ballistic nomenclature does not include the a, 
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components of force and moment. However, the relation "between the aero- 
dynamic factors and the constants of the |(x+) equation given by equa- 
tion (132) shows clearly that should be added to mq in the inte- 

grated equations of motion as they apply to flight in the aerodynamics 
range. !^rthermore, since + Zq) is the force corresponding to the 

moment (m^ + mq) just as Jg is the force corresponding to the moment 
Jj[, should be added to Zq in order to have the equation for Jg 

consistent with the equation for Jg. 

The lift and drag factors are introduced in the same manner in the 
ballistic development as in this report (see eq. (l9))^ that is 


hence 



^oL " 



L 


( 1 ^ 0 ) 


Since they may be of interest in reading the ballistic literature, 
the equations between the aerodynamic and ballistic coefficients are 
listed below: 


% = 






K, 


■N 


S 

2d^ 



Ks 


i i 
F d^ d 




- If (d) 


Kp 


1 ^ I p 
“ f d^ ^Pp 


bfT - - 1 i (a) 


KxF 



%r 





and 


Kl = 


s 

2 d^ 



(i42) 


The x"''(t’'') history is given by equation (j^); however, a more suit- 
able form for the purposes of comparison is given during the development 
of equation ( 7 ^) Appendix A, as follows:' 


m 


dV 

dt 



(l!+ 3 ) 
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Equation ( l43 ) becomes , 
tion ( 68) , 


after transforming from time to distance by equa- 


The comparable equation from reference 20 is given on page 8 of the 
reference as follows: 


dz 



(145) 


where 


V = velocity of projectile (= V) 
z = distance along range ( = x) 



If equation (li+if) is transformed into ballistic notation^ it becomes, 
using Cp = as specified by equation (l9)> which is correct to first 

order since = 0 for projectiles and missiles, 


dV V ^ 

d^ = ■ pT D 


(11^) 


Comparison of equation (lk6) with equation (ik^) shows that the results 
of this paper agree precisely with those of reference 20. 


The cp(x'*') history is given by equation (86). The comparable equa- 
tion from reference I 9 is equation (l^J-) of that report: 


qp" + CjCp* - Cg = 0 

where (in the notation of ref, I 9 ) 


K- 


I 


% 

m 


^2 


(ihT) 


(148) 


A study of the nomenclature of reference I 9 discloses the following rela- 
tions between its nomenclature and that of the present report: 
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9 = cp 



(149) 


After it Is transformed into the nomenclature of reference I 9 by the use 
of equations (l48) and (l4-9)j equation (86) becomes 

cp" + Cicp‘ - C2 = 0 (150) 

Comparison of equation ( 15 O) with equation (l4?) shows that the results 
of this paper agree precisely with those of reference I 9 . 

The differential equation for the A(x''') history is equation (93) > 
whose constants are given by equations (99) through (l02). The comparable 
equations from reference 21 are equations (4o) and ( 5 ) of that report: 

S" = (Jl - ivjp) I + ±{vj^ + lJg)^« - ^ (151) 

^ + l^y = k^ei'Pi + k (152) 

where k^_, k^, and cpg are given by equations (6) and ( 7 ) of refer- 
ence 21: 





tjrr 

Jl - ^ - (2Jt 


Jl + j 




(Jp - md^J^/A) ' 




1 

d 


dz 


(153a) 
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Experience indicates that the term i(vJxp + iJg)^' in equation (151) 
is negligible compared to the remaining terms ^ and consequently the bal- 
listic equation for the transverse displacement reduces to 

S" = (Jl - ivjj.) 1 - if (155) 

d 

A study of reference 21 shows that its nomenclature is transformed 
into the nomenclature of this report by equations (13T)> (l39)j s.nd (l^i-O) 
together with the following: 

'1 

S = - pZA. 

C = I * (156) 

Z = X 

- 

If equation (8b) is transformed into the ballistic notation^ it 
becomes 

s" = (J, - ivj,) 1 . i| . sisiniJio) 

where Uq = u at t = 0. 

Comparison of equation (157) with equation (155) shows that the 
transformed A(x"^) equation agrees exactly with the ballistic equation 
except for the term 2ig(u - Uo)/u®. For representative tests in the 
aerodynamics range, this term is negligible compared to .the remaining 
terms in equation (157) * Hence, it may be concluded that the agreement 
between equations (157) and (155) is satisfactory. 
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From the comparison of equation. (93) for the |(x'^) history with the 
corresponding ballistic equation (152), it can be seen that both equations 
have the same form. Hence ^ the question of agreement between the two is 
concerned with a comparison of the constants of the equations. The aero- 
dynamic factors are determined solely from the constants <3^ and d4 (see 
eqs. (93) and (l32) through (l35)) and it is believed sufficient for the 
purposes of this report to confine the comparison to these two constants . 


Comparison of equation (93) with equation (l52) shows that the real 
and Imaginary parts of the exponents in equation (93) are related to the 
exponents in equation (152) as follows: 




X 


+ 


2d 



Jl + Jp 



md^JjY 

A 


+ 


1 

a 



(2Jt - Ja) 


md^ 

A 


dz 


(158a) _ 


d4R x+ 



Jp + - 


k2 



1 

a 




(2Jji - Ja) 


md^ 


dz 


(158b) 


^2I x+ = cpio + 




vA(l+g) 

B 


dz 




dz 




( (159) 


j 


Differentiating equations (158) and (l59) 'with respect to z, noting from 
equations (156} that x = z, and limiting the comparison to t = 0 gives 


2d 

pi 


'‘2R 



- J. 


D 




md^JA 

A 


+ 


1 

a 


Jl + Jd - p 




_ J 


(l60a) 
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2d 


^ 4 R 


Jl 


Jd 


+ 


*% 

k 2 


+ 


1 

02 




1 

0 


‘^D 




(1601)) 


v-l 


^ 2 I 


= ^ (1 
B ^ 


o) 


> 


\xl 



(1 


a) 


(161) 


Up to this point explicit expressions have not been derived for the 
real and imaginary parts of dg and d4_, since such expressions are not 
required for the analysis of data from the range. However, they may he 
readily obtained from equations (132) through (135)* Ths solution given 
below is valid for all values of s or 0 except the region for which 
0 < s +1 for which 0 is imaginary. If mo, is negative, correspond- 

ing to 8 < 0 or 0 > 1, the model will have arrow stability. If m^;, is 
positive but s > 1, corresponding to 0 < 0 < 1, the model will be gyros- 
copically stable. Hence it is believed that the following solution will 
cover all cases of practical importance 



vK(l+o) 
•ial - 2 

vK(l- 0 ) 
<i4l = 2 


(163) 


Transforming equations (162) and (163) into ballistic notation and 
multiplying by 2 d/(j.Z, as required for comparison with equations (160) 
and ( 161) , gives 
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2d 

M-l 


dgR - - (J^ - Jx) + i + Jd - k"^JH “ 2 


md^ 

“T '^T 


(l64a) 


— <3.4R = - (Jl - Jd + - i Tjl + Jp - k'^J^ - 2 S|_ j,j 


(l6kb) 


2d . 

7l 


2d . 

ITT ^4i 


f (1+a) 


vA 

B 


(1-d) 


(165) 


Comparison of equations (l64) with equations (l60) and equations (l65) 
with equations (l6l) shows that corresponding equations agree except for 


the term 



md%^ 

■“a 


which is missing entirely from equations (l64). 


Differences in the derivation of the equations of motion make it difficult 
to ascertain the reason for the discrepancy. ~A careful examination of 
the terms neglected in this paper in linearizing the equations of motion 
suggests that the term in question is actually a second-order quantity. 

In practice , the extra term is so small for the majority of tests carried 
out in the range that equations (l6U) can he said to agree satisfactorily 
with equations (l 60 ). 


To summarize, the results of this paper are in good agreement with 
the results of ballistic theory. The equations of motion are the same 
except for an additional term in the equations giving the transverse dis- 
placement and for two additional terms containing a common factor in the 
relations between the aerodynamic coefficients and the damping rates of 
the oscillatory components of the angular motion. It is believed that 
these additional terms are negligibly small \uider the majority- of repre- 
sentative test conditions in the aerodynamics range. 


CONCLUSIONS 


Equations of motion have been developed for conventional aircraft on 
the one hand, and for rotatlonally symmetric missiles and projectiles on 
the other. Similar mathematical derivations and standard NACA nomencla- 
ture were used for both developments. The essential difference between 
the two cases lies in allowing the axial component of the angular velocity 
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to be large In the projectile and missile case. Otherwise^ the conditions 
of flight are limited to small changes in velocity and orientation^ as is 
customary in treating the first-order dynamic stability of aircraft. 

Two novel features appear in the derivation of the equations of 
motion. First, the aerodynamic coefficients are introducted as a formal 
series expansion in the components of the linear and angular velocities. 

This approach allows the Introduction of coefficients of quadratic terms 
involving the axial component of angular velocity, which are second-order 
for the aircraft case but first-order for the missile and projectile case. 
Second, it is shown that the aerodynamic forces apd moments are independ- 
ent of the orientation in roll for rotationally symmetric missiles and 
projectiles. Advantage is taken of the Independence to place the coor- 
dinate system used in the missile case with one axis along the missile’s 
axis of rotational symmetry and another axis in a plane fixed in space. 
Hence, the position of the coordinate system is independent of the roll 
orientation of the missile. 

Criteria are derived for the dynamic stability of missiles with arrow 
stability and projectiles with gyroscopic stability. It is shown that 
spin-stabilized projectiles may be completely dynamically stable, a result 
in agreement with experiment. On the other hand, it is clear that the 
dynamic stability of spin-stabilized projectiles is a delicate balance 
between lift, damping, and Magnus forces and that all three elements are 
required to predict adequately the character of the motion. 

The equations of motion are applied to the flight testing of rota- 
tionally symmetric missiles and projectiles in an aerodynamics range. It 
is shown that the flight records obtained in this facility can be analyzed 
to give all the aerodynamic coefficients required to predict the first- 
order dynamic stability of the missile or projectile. 

The results of this paper are compared with the accepted aeronautical 
theory for aircraft and the corresponding ballistic theory for missiles 
and projectiles. The agreement is considered to be satisfactory throughout. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 1, 
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APPENDIX A 

LINEARIZATION OF THE EQUATIONS OF MOTION 


In the author’s opinion^ the proper procedure to be followed in 
linearizing each of the equations of motion is to separate the exact equa- 
tion into first- and second-order terms and to place all of the first- 
order terms that are to be retained on the left-hand side and all of the 
second-order terms that are to be neglected on the right-hand side. This 
procedure allows a term-by-term comparison of the second-order terms 
neglected with the first-order terms retained at one step in the deriva- 
tion. As a result, it places the neglect of second-order terms on a firm 
basis that can be re-examined readily at any time should there be some 
doqbt as to the relative magnitude of any particular term. 

The procedure just recommended is carried through only in part. All 
of the second-order and some of the third-order inertia terms are included. 
However, only part of the second-order aerodynamic force and moment terms 
are written down. Although it would be desirable for the sake of rigor 
to include them, their omission is justified in order to make the deriva- 
tion more concise since there is rarely any doubt as to their magnitude 
relative to the first-order aerodynamic terms. It is the comparison of 
the Inertia terms with the aerodynamic terms that is ordinarily subject 
to scrutiny. 

No explicit limit has been set so far on the values of the aerody- 
namic coefficients themselves; yet, there is the implicit assumption that 
their magnitude ranges in value from unity to ten. This restriction can 
usually be met by the appropriate selection of representative length and 
area. Also, there is the implicit assumption in the missile and projec- 
tile case that the axial component of the angular velocity (the spin) be 
bounded. Again, this restriction is normally satisfied by the values of 
the angular velocity which actually occur in practice. 

It will be noted that an anomaly exists in setting limits on both 
^ and p,q,r for a given aircraft. Once the physical and aerodynamic 
properties of the aircraft are specified together with the initial condi- 
tions of flight, the values of p, q, and r are determined. Physically, 
the consequence of assigning limits to both cp,0,t and p,q,r is to assign 
a certain range of values to the physical properties of the aircraft, that 
is, to its moments of inertia, weight, and center of gravity. As a result, 
it is desirable to examine the values of p, q, and r occurring during the 
flight of any particular missile or projectile in order to make certain 
that their magnitudes fall within the limits of the theory. 

In the following sections, the equations of motion will be developed 
up to the point of listing the first-order terms retained on the left-hand 
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Bide and the second-order terms neglected on the right-hand side; but a 
numerical comparison of their relative magnitudes will be omitted. In 
the aircraft case, the terms neglected are those customarily neglected. 
In the missile and projectile case; numerical examples were worked out 
for representative missiles and projectiles and it is believed that the 
terms neglected are truly second-order for the majority of conditions 
that will occur in actual practice. 


Aircraft Equations With Time as Independent Variable 


X component of force . - This equation is derived by resolving the 
vector force equation (3) aJ-ong the space-fixed Xq axis . 

Since 


M= mV (A1) 

the component of equation ( 3 ) along Xq becomes 


m cos'^ (XXo) + cos^ (YXq) + 

^AZ <=os-^(ZXo) " “€ 7o (A2) 

Substituting equations (ll) and ( 15 ) for F^^Y;Z equations (l) for 

cos ^(XXq); cos -^(iXo)? and cos -^(ZXo) in equation (a 2); and at the same 
time neglecting all second-order terms on the ri^t-hand side of equa- 
tion (A2) ; gives 


dV ^'^o 

mfcos ^VXo - V sin^VXo 


dt 



CXo 





CX 


31. _ 

q 2V 


Cy 0 


- mg sin 7 q 


(A3) 


cos'):VXo being approximated by 


cos '^^VXq = 1 


slng^VXo 

2 


and the velocity; V, is given by 
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V = Vo + u 


> 



(A^) 


If equations (a 4) are expanded In a series^ equation (A3) may be separated 
into first- and second-order terms with the first-order listed on the 
left-hand side and the second-order on the right-hand side, as follows: 


du p7o S 
“ dt ^ 2 


+ Cy d + Cy . oTT — + Cv — + C Y 6 +2 tt~ Cy j + 
-^o a 2 Vq -'^q. 2 Vq Vq ^o/ 


mg Bin j-o = I ff’ ulr? -J: + mV ain^VX^ 


r„ pVo^lS / 2u ^ 

- — ^ * 


al u 


ql 


2u 


it at vt - ^ 


pvfs 


Hi 
2V 


(A5) 


If the right-hand side is neglected, the left-hand side is divided by 
(pSVo^), and the equilibrium conditions (eq. (l8)) and lift and drag 
coefficients (eq. (19)) are introduced, equation (A5) becomes equa- 
tion ( 20) . 

A slight variation of equation (A5), which will be useful in deriving 
other equations of motion, comes from equation (A3) hy introducing the 
equilibrium condition of equation (18) 

^o = - 


namely. 


m 


dV _ pSV2 

dt 2 




ai . 


a 2V 


2V 


CZo® 


(A6) 


Y component of force . - This equation is derived by resolving the 
vector force equation (3) along the Y axis. 

_The components of M and V are given by equation (6), the components 
of are given by equations (ll) and (15), the components of Fq are 

given by equation (l?)^ and the transformation of moving axes is given by 
equation (9)- Equation (3)^ modified by the foregoing, may be resolved 
along the Y axis as follows: 
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m 


dt 


sin 3 + V ~ cos 3 + rV vl - sin^a - sin^P - pV sin a 


pV^S 


Cy^P + Cy. H + Cyp 11 + Cy^ + mg(t sin 7o + 9 cos 7 q) (^7) 


Tile following quantities are approximated by 


cos p = 1 ^ 


sin 3 = p 




sin^a - sin^P = 1 - ^ — i-^! — I 

2 


The values of r and p are given by equation ( 7 ) 

dt 

r = — - cpq 

P = - 0r 

Q*b 

After division by V and the use of equations (Ab) and (a 6) for V and 
dV/dt, equation (A7) may be separated into first- and second-order terms 
with the first-order listed on the left-hand side and the second-order 
on the right-hand side, as follows: 


m 


I - i) - ^ (V%- ^ 


+ 


r(a^+B^) 


9 cos 7 q) - p + 9 q + — ^ ' + paj + 3 f Cy^a + Cy^ ^ + 


dZ 


qZ 


_ £§V 


CXq ^ ^ I^Cype ^ + Cy^cp + 


Tl 


ql 


u 

Vo 


pSVo 

CypP 


ing / ^ 

— (,t sin 7o + 9 cos Jq) 


(a8) 
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If the right-hand side is neglected, the left-hand side is divided by 
(pVqS), and the lift and drag coefficients (eq. I 9 ) are introduced, equa- 
tion (a 8) becomes equation (2l). 

Z component of force .- This equation is derived by resolving the 
vector force equation (3) along the Z axis . 

The components of TK, V, F^, and Fq are given by equations (6), (ll), 
(15), and (17), respectively. After transformation to moving axes by 
equation (9) and division by V, equation (3) resolved along the Z axis 
becomes 


/ tia 

ml ^ cos a + 


dV sin a 
dt ~ 


+ p sin p-q - sin%, - sin^p 


) 


isy 

2 


-CZo - CZa,a - ^ il) ^ (- ® ’'o + 7o^ 


(A9) 


The following quantities are approximated by 

sin a = a 




1 - sin^c6 - sin^P = 1 



cos a = 1 


or 

•F 


The value of q is given, by equation (7b), namely 


q = — -t-qjT 


By the use of equations (a 4) and (a 6) for V and dV/dt, equation (A9) may 
be separated into first- and second-order terms as follows: 
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m 




m 

Vo 


u 


- e sin 7o + ( 1 - ^ ) cos /o 


= m 


— + cpr - pP - ^ — 2 ^ 


pVoSa 


0,1 ql 


J + “2 [^^o^ "■ ^X^2V ^ ‘^Xq 2V ■*■ V 


i(. a + e Sin ro) - (cz,T H ) 


(AlO) 


If the right-hand side is neglected, the left-hand side divided hy (pVqS), 
and the equilihrium. conditions (eq. (18)) and lift and drag coefficients 
(eq. (19)) are introduced, equation (AlO) hecomes equation (22), 

X component of moment .- This equation is derived hy resolving the 
vector monient equation (4) along the X axis « 

The components of H, and M are given hy equations (6), (T)^ 

(12), and (15). After transformation to moving axes hy equation (9), 
equation (li.) resolved along the X axis hecomes 


^x ^ ' "^xz ^ - ly) - qp J; 


xz 


H + C, C7 


- ^ 2V ^ ^ 2V 


(All) 


The values of p and r are given hy equation (?) ^ namely, 

n - fir 

dt 
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The value of V is given by equation (a 4)^ namely^ 


V ^ u 

Vo Vo 

Equation (All) may be separated into first- and second-order terms as 
follows , 

d^cp pSZVo^ (^P p , p 2 , p, . .1 ^ ^ n, dt 

% ^ - JXZ ^ ^ ^ dt ^ 27o dt + ^^r 2Vo dt, 


g 

= Ix (r0 + 6r) - Jxz ( q<p + Cp4) - qr(l 2 - ly) + qp Jxz + ~ ^ + 


Cl ^ + C 2 

3 2Vq ^P 2Vq ^r 2Vo 


C, 0 £A + C7 e ^ 
‘•p 2V ^r 2V 


(A12) 


If the terms on the right-hand side are neglected and the left-hand side 
is multiplied by p/%^ pSZVq^, equation (A12) becomes equation (23). 

Y component of moment .- This equation is derived by resolving the 
vector moment equation (4 ) along the Y axis . 

The components of H, w, and M are given by equations (6), il) > 

( 12 ), and ( 15 ). The value of q is given by equation ( 7 ) to obtain q, 
namely 

a = 

^ dt= ht 

After transformation to moving axes by equation ( 9)7 equation (4) resolved 
along the Y axis becomes 

ly + ^p(Tx - Tx,) + Jxz (P^ - ^2) 

= ^^2^ ^ + ^mq (A13) 


The value of V is given by equation (a 4), namely, 


V u 

TT- = 1 + 7 ;— 
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Eq_uation (AI 3 ) may be separated into first- and second-order terms as 
f ollo-ws : 


lY 


d^e pszvo^ 


dt2 


2 y '^mopt- + ^d, 2 Vq 2V(; 


= -lY + rp(lz - Ix) + Jxz(r2 - p2) + ^ (^2Cnio + + 


Cm, 


0,1 

^ £Vo 


+ Cm, 


q. 2Yr 


-Si- 1 + PC^V 


Cmq^9 || 


(Al^) 


If 'the right-hand side is neglected, the left-hand side multiplied by 
u/pSIVq^Ky^j and Cm^ set equal to zero in accordance with equilibrium 
conditions specified by equation ( 18 ), equation (A14) becomes equation {2h) . 

Z component of moment .- This equation is derived by resolving the 
vector moment equation ( 5 ) along the Z axis. 

The components of H, co, and M are given by equations (6), (T)!? 

( 12 ), and ( 15 )* The values of p and r are given by equation ( 7 ) to 
obtain p and r, namely, 

• _ d ^9 _ d(6r) 

^ dt^ " dt 

. d^ijf d(cpq) 

“ dt2 " <at 

After transformation to moving axes by equation (9)^ equation ( 5 ) resolved 
along the Z axis becomes 


^Z 



d(cpq) ■ 


d^cp 

d(0r) ■ 

_dt^ 

dt 

- Jxz 

_dt^ 

dt ■ 


+ Plly " iCp^X " ^Jxz) 


pSZY‘ 


CnpP + Cn^ ^ + C^p 2 ^ + Cn^ 2 ^ ) (AI5) 


rl 


V is given by equation (a 4) , namely, 


V, 


V ^ u 

= 1 + rr- 


The value of 
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equation (AI5) may be separated into first- and second-order terms as 
follows , 


^ d^ ^ d^ pSZV'o'^ f ^ a . M , n J_ ^ ^ n„ ^ 

iz ^ - Jxz ^ ^ H 2Vo dt + 2Vo tt ^ 2 Vo at 




a(6r) 


Iz - Jxz + P<lUx - Iy) - Jxzgr + 


pSlVo-" U 




pz 


+ Cr 


pZ 


2Vo ""P 2Vc 


rl ^ pS^Vo 


rl 


ql 


2Vq / “ 2 V “p® 2Vo ’*' 2V, 


:) 


(A16) 


If the right-hand side is neglected and the left-hand side is multiplied 
by n/pS2Vo^K2^^ equation (A16) becomes equation (25). 


Symmetrical Missile and Projectile Equations With 
Time as Independent Variable 


The missile and projectile equations are derived by resolving the 
vector force and moment equations (3) and (4) along the XYZ axes, except 
for the X component of the vector force equation, which will be resolved 
along the space-fixed axis Xq instead of the body axis X. The deriva- 
tion -VTlll proceed as follows: First, the component of the exact vector 

equation will be written down transformed to the moving XYZ axes, except, 
of course, for the Xq component equation; second, the kinematic rela- 
tions will be substituted for certain of the dependent variables , various 
functions will be expanded in series, and the resulting equation will be 
written down separated into first- and second-order terms with the first- 
order terms on the left-hand side and the second-order terms on the right- 
hand side. 

The following relations will be used: The transformation from space- 

fixed to moving axes is given by equation (40b); the direction cosines of 
the angles between Xq and X,Y,Z are given by equation (l); the components 
of M and V are given by equation (6) and of H and to by equation (39); 
the components of are glv£n by equations (ll), (29), and (42), of 

Fq by equation (44-), and of M by equations (l2), (29), and (42); the 
magnitude of the velocity, V, is given by condition (2) preceding equa- 
tion (2) and the axial component of the angular velocity, p, by equation 

( 34 ). 
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The following series expansions will be utilized: 

sin P = P 


Bxn a = a 


cos a = 1 - 


g,^ ft T P^ 

— COB p = 1 - 

2 2 
2 


■Jl - sin^ a - sin^P = 1 - 


+ P 


wT“ 1 - 

^oy Vo 


(AIT) 


In deriving the Y and Z components of the force equation^ the 
following approximation will be used 


1 OT 
V dt 


% ^ 


(AI8) 


As will be shown in the derivation of the X component of the force equa- 
tion, the above equation neglects the terms pSV/ 2 m and g sin 7 oA • 

From the limited experimental evidence available, it. is believed that the 
Cx term is truly negligible under all practical circumstances . On the 

other hand, the gravity term is strictly negligible only for horizontal 
trajectories or in those cases for which the drag force greatly exceeds 
the weight. In other cases, for example, a bomb falling along a steep 
trajectory, the gravity term should be included. The gravity term will 
be omitted from the present treatment since the application covered in 
this ■ paper is to the aerodynamics range in which the trajectories are 
nearly horizontal and the drag force is many times the weight in the great 
majority of tests. 

X component of force. - 


dVv 


m 


dt 


= Fj^cos + Fj^jCOB ^XqY + Fj^cos ^XqZ - mg sin 


(AI9) 


where is the component of V along Xj^ and is given by 


Vx^ = V cos^-^ 


m 


(cos^VXo) (Sj “ (V sin^VXo) 

_ pSV^ 

“ 2 


'd^ VX( 
dt 


- CXo - CXpp (^|i 


- mg sin 7 q (A 20 ) 
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du . 



2 ^X 0 ^ - 2Cxpp 






mg sin 7o 


= m 


sin2 ^VXr 


f)MV Bln^VXo)^'"^° 


pVo^S 




u 


+ Cx- 


2II 

PP V 2Vn 


- ^Xpp^ ~ 


(A21) 


Dividing equation (A21) by (pSVo^) and neglecting the right-hand side 
gives equation (46). 

A variation of equation (A21), which will be needed in the deriva- 
tion of the Y and Z force component equations, derives from equa- 
tion (A20). In this variation the Cxpp and g teims are neglected; 

hence equation (A20) may be written 


m 


5Y 

dt 


+ Cx, 


pSV2 

2 


1 ^3 In^ ^ 



+ (V sln-):VXo) 



(A22) 


If the right-hand side is neglected, equation (A22) becomes equation (A18) . 
after the terms are rearranged and the equation is divided by mV. 


Y component of force.- 


dMy 

— -1- rMx = Fay + ^ 


GY 


(A23) 


m 




psVo T/ 

Vdt 

dty 

2 1_\ 


CYu^-OxJP + CYj II + Cy^ 5U |i + M „ + 


Po^ I da 


= m 


dt 


Po^ 

1 

de 

2Vo 

2Vo 

dt 

-t- r 


ras 

2 


2Vo dt "^r 2Vo dt " 2Vo 

^mg Bin'7Q 
■ \ ^ 


pSVo u 




CYp + CXo)P - 


2L _J_ iia r ±L JSJl 
"^X^p 2Vo 2Vo dt " "'Xqp gVo 2Vo J 


pSVp p«2 


2- -^ 


efe if + “X 


_iL 

qp 2 Vq 


mg sin 7 qV,, u 


Yr 


Vo 


(A24) 
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After division by (pSVq) and neglect of the right-hand side;, equation (A2ij-) 
becomes equation (47). 

Z component of force . - 


dM7 

"dt '1^X= FaZ + Fq2 (A25) 




aak - hZ. u O ( n o 1 ft n ^ , n i ciy 

<34 dty 2 _ 2 Vq dt 2 Vq dt 


I d0 


2Vo 


Po^ Q n T' 

P + Cz, 


^PP 2Vq 2Vq dt 


+ Cz M _J_ M 

2Vq 2Vq dt 


- ^ (cos 7 o - 0 sin 7n) 


= m 


f da 
dt 


qp^ 


pSV, 


o u 


2 Vr 


-Cz + Cx_ ] g + 


G — _L ^ ^ P rZ 

2pp 2Vo 2Vo dt ^ "^Zrp 2Y^ 


pSVp P'Z 

2 2Vo 


n ft 4 - p ^ , 

Czpp 3 Cz^^ 27^ at 




^2rp 2Vo y 7o Vo ~ 


(A26) 


(pSVq) and the right-hand side is neglected, equa- 
tion (A26) becomes equation (48). 

X component of moment . - 


dHx 

— + qHz - rSy = % (A2T) 


Since rotational symmetry requires that 
becomes 


Iz = lyj equation (A2T) 


d2(p pSZVo" 




I 

h 2Vc 


dcp 

dt 


Pq3 


U 


2Vo Vn 


pSZVp^ / p» X 

2 V Zp Vo 2Vo 


C 7 0 
'■P 



+ 


d(0r ) 


dt 


(A28) 


If it is multiplied by p/pSZVo%x^ the right-hand side is neglected, 

equation (A28) becomes equation (49)* 
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Y component of moment. - 


dHy 

dt 


+ rHy = My 


(A29) 


d^e T dt ^ + 


^ dt ^ 2Yo dt 


■'nil 


PP 


Po^ R ^ n P^ _L ^ 4 - P Pq^ ^ 

2Vo ^ 2Vo 2Vo dt “rp gy^ 2 Vq 


- lyrp * + ^ (^2CmoA' + Cm^ gy^ + 0^^ 2^^ 


aZ 


ql 


n \ , pSZVp^ p«Z 

Sp^V 2 2Vo 


Z ^ , p r Z 
CmppP + 2Vo dt ^ '^“rp 2Vo 


^ (A30) 


Through multiplication by p/pSZVQ^Ky^ and neglect of the right-hand 
side, equation (A30) becomes equation (50). 


Z component of moment. - 


dH^ 

^ - q% . Mz 


(A31) 


iJ - ’^XPO at 


d6 pSZVo 


CnoP + ^n. ov~ 41+ 


Z df 


2 " '"’^p 2Vo dt 2Vo dt 


Po'^ 


Po^ Z da 


Po^ Z d0 


^ 2Vo 2Vo dt ^ ^^qP 2Vo dt 


= lyqp’ + ^ (2CnpP + ^n^ ^ + Cnr ^ + C 


^ap ^Vq 


a + 


pS^Vo^p'Z „,p JL 

2 WZ 2Vo dt ^ ^=^qP 2Vc 


(A32) 
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Through multiplication hy p/pSIVq^K 2 ^ and neglect of the right-hand 
side, equation (A 32 ) becomes equation ( 51 ) • 


Symmetrical Missile and Projectile Equations With 
Distance as Independent Variable 


The missile an& projectile equations with distance as Independent 
variable are derived in a similar manner to those with time as independent 
variable. The procedure followed closely parallels the preceding section. 
First, the component of the exact vector equation will be written do-vm 
and for those components which are resolved along the XYZ axes it will 
be written transformed to moving axes. Second, the kinematic relations 
will be substituted for certain of the dependent variables , various fimc- 
tions -i-rill be expanded in series, differential quantities will be trans- 
formed from time to distance as Independent variable, and the resulting 
equation will be written down separated into first- and second-order terms 
with the first-order terms on the left-hand side and the second-order 
temis on the right-hand side. 

The same relations will be used as in the previous section. In addi- 
tion, the transformation from time to distance is given by equation ( 70 ). 
Also, since the principal application of the equations with distance as 
independent variable is to the aerodynamics range, the earth axes will be 
oriented as they hormally are in the range with Zq vertical and Xq 
horizontal j consequently, 7q = 0. 

It is evident from the relation, V = dx/dt (eq. (6?)) aiid the develop- 
ment of equation (A22) that the X-component force equation gives the 
relationship between time and distance. Consequently, it is proper that 
time should remain the independent variable in this equation. However, 
the form of the x'*‘(f*') equation (T^) may not be familiar and the steps 
leading from equation (A22) to equation (7^) will be presented in this 
section. It should be noted that in this development the term is 

neglected . 


Application of this theory to the aerodynamics range requires equa- 
tions for the transverse displacement of the trajectory, y and z. The 
equations for y and z are derived by resolving the vector force equa- 
tion ( 3 ) along the space-fixed Yq and Zq axes. Their development will 
be presented in this section. It should be noted that the only forces 
retained in deriving the y, z equations are lift, drag, and Magnus 
forces (CL^,CDo.Czpp)- 

Time-distance equation .- If equation (A22) is multiplied by dt/mV^ 
and the right-hand side is neglected, it becomes 


dV 


2m 


dt 


(A33) 
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Since V=YQ&tt=0, the integral of equation (A33) may be written 


V = Vr 


1 + 

2m 


(A3^^) 


If equation (a 3^) is substituted into equation (A33) and the relation 
V = dx/dt is used, the resulting equation may be written 


dt^ 


+ 


psvp^ 

2m 


1 + C- 


pSVc 

’^o 2m 


-2 


= 0 


(A35) 


After multiplication by m/pSVo^, equation (A35) becomes equation (74-). 
Transverse displacement of trajectory (y, z equations) . - 


l-I = Fxcos ^YqX + Fycos + Fzcos 

dt^ 

d^z 
dt‘ 


m 


m = Fxcos + Fycos + Fzcos 


(A36) 


m 


m 


^ = psy^ f - C llr - R + C-7 r. 

T!a ? V + ^Zpp 2Y a- 


dt 


d^z pSV^ n r. n 2^ 0^. 

= — Cv S - C .7 a - Cy- ^ p ) + mg 

^p2 2 V o ^cf, ^pp 2v ^ J 


(A37) 


Transforming equations (A37) from time to distance, dividing by V^, 
and separating into first- and second-order terms gives 


m 


^ f CzJ - C2, 


dx" 


2 dx 


a 


•■PP 2Vr 


a + Cx t 




— U' u. I oxr n\T ir 


V2To ev„ Vo 


,2 2 dx 2 V Za Zi 

dx ^ 


m 


Pq^ 


"PP 2Vq j " Vr ^ 


mg 


) (a38) 




PP r /£li. _ ^ Ji. ^ - 2 .a 


2Vo 2 Vq Vo y 


V, 


2 Y. 


O w 
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For small angles 




dz 

dx 


a 


e 


if lift and drag coefficients are introduced and the right-hand sides are 
neglected; equations (A 38 ) become 


, ps 

dx^ 2 


m 


d^z 



- %P 

(cl„«.+ Czpj,!^ p) 


= 0 





(A39) 


After multiplication by (pS)“^; equations (A39) become equations (75) and 

(T6). 


Y component of force . - 



+ rMx - ^AY + ^GY 


(A23) 


m 



+ 






1 dp 

2 ^ ^ 


^Yr 


1 ^ 

2 dx 




a + 


Poi i da , p I de 

% 2Vo 2 dx ^Yqp 2 V 0 2 dx . 


= m 


2 \dx dx / 2 dx 


pS u 
2 Vo 



_pL 

Yap 2Vo 


a-f- 


Cy. 

^ap 2Vo 


I ^ , P qZ ^ 

2V dt Yqp 2 Vq 2V J 


£§ 

2 


pU 

5vY 




Z da ^ 

2V dt ^Y, 


iL 

qp 2V 


(a 4 o) 


Dividing equation (a'+O) by (pS) and neglecting the right-hand side gives 
equation (77)* 
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Z component of force . - 


dMz 

dt 


- q.Mx = FaZ + ^GZ 


(A25) 




dx 


'PP 2Vc 


Po^ ^ d£ „ Po^ I ^ 
__ P + O7,. o-ir o Zrp pv^ 2 dx 


mg 

^ 2 dx 2Vo 2 dx J Vq^ 


pp 


m 

2 


(- to _ ae 


dx dx 


dx 


pS u ±i±ip, 

2 Vo V 2Vo 2 Vo 2V dt ^ 


yl rl\ _ £§. £jli. f Q„ 3 + C2. ~ + C7 ^ (a 4 i) 

'^Zrp J 2 2Vo ^ ^PP^ ^3p 2V dt Zrp 2V J Vq^ Vq 


Dividing equation (a4i) by (pS) and neglecting the right-hand side gives 
equation (iB) . 

X component of mment . - 

^ + qH^ - rHy = % (A2T) 

Q. Xf 

With the introduction of a constant rolling moment, equation (A27) 

becomes 


ix 



p3Xp ^ 

m dx 



+ C 


I 

^p 2 dx/ 


Ix d(er) psi rl 

y2 dt 2 2V 


(Al^2) 


Multiplying equation (a 42) by p/pSlKy^ and neglecting the right-hand 
side gives equation (79)* 

Y component of moment. - 


_ + rHx = My 


(A29) 
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^ d^e . ^ Po dt psz r 2 <to , I de 

TT + ^x ^ di ■ 2 ^ 2 35 E i ^ ^ d 5 E 


Sp io ^ ■" '^B 


^ ^ _2 ^ 
pp 2Vo 2 d^ ^ 2Vo 2 dx 


_ ^ ^ /p^ p»Z\ pSZ _p^ JL Rj-P _L^4. 

2 dx \7o Vo VoJ 2 2 Vo Vo % 2V dt 


rz\ pSZ p'Z 


^rp 2vJ 2 2Vo V^Pp ^ 


Ju ^ 

2 V dt 


p El 
2 V 


(A^ 3 ) 


If equation (a^ 3 ) is multiplied by fi/pSZKy^ and the right-hand side is 
neglected, it becomes equation (80) . 


Z component of moment. - 


q% “ 


d^V ^ Pq pSZ ^ ^ 2 dp , 

_ - Ix ^ ^ — I^CnpP + Cn. 2 to + 


Pq 2 2 da 


Pq 2 2 d0 


2 dx ^op 2 Vq “^“ap 2 Vq 2 dx + 2Vo 2 dx 

= _I^.^^l_ + ^X£ 2 £ll_e^£ 2 _ 3 ^An a + 

X Vo dx Vo V Vo 2 2 Vo Vq 


(A 31 ) 


P _L 1 ^+P £l^4.££llll 

2V dx ^qp 2 YJ 2 2V^ 


(^ap“- + ^^qp H) 


If equation (A 2 t- 4 -) is multiplied by p/pSZKy^^ and the right-hand side is 
neglected, it becomes equation (81). 
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APPENDIX B 


MODIFICATION OF THE SYMMETRICAL MISSILE AND PROJECTILE EQUATIONS TO 


INCLUDE THE EFFECTS OF SMALL AERODYNAMIC ASYMMETRIES 


Forces and moments due to aerodynamic asymmetries are not independent 
of the roll angle as are the principal aerodynamic f£rc£S and moments. 
Hence, they -will be defined as components along the 1 , 2 , 3 body-fixed 
axes as follows: 


1,2,3 


M 


1,2,3 


‘^Fi, 2,3 



(Bl) 


The components of Cp and Cf^ along the 
CXq (due to asymmetries) = Cpi 

CYq = CpaCOB 
= CF3CO8 

^lo = 

Cjjj^ = Cj^gcos 
= Cj^COB 


X,Y,Z axes are given by 


cp - CpgSin cp 
(p + CpgSin (p 

> 

cp - CMgSin cp 
cp + CMssln cp 


(B 2 ) 


It is assumed that the coefficients, Cp3_^2,3 *^Mi,a,3 first- 

order quantities, in contrast to all other aerodynamic coefficients. This 
is evidently the case if the asymmetries are small. Also, if Cp or C14 
are due to control-surface deflections, 5 , they will have the form Cp^S 
and, hence, the requirement that they be small is equivalent to assuming 
that the control surface deflections be small. Consequently the variable 
velocity, V^, may be replaced by the constant, Vq^, since the term 
2 C u/Vq will be second-order and hence negligible insofar as the line- 
arized theory is concerned. 


Certain of the models tested in the aerodynamics range have small 
asymmetries of the type Just described. It is desirable, therefore, to 
derive the modifications to the equations of motion which are required to 
account for the effects of the asymmetries Two of the components have 
already been accounted for. The X force component, Fi, is a constant 
and its coefficient, Cp^, may be considered to be Included in the axial 
drag coefficient, CXq. The X moment component, M^, has been included in 
the 9 (x'^) equation by the Iq term. It should be noted that in this 
particular equation the Cmi coefficient is not required to be small 
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(first-order). The remaining components affect the a(x'^) and |(x+) 
motions^ and the changes to the equations involved -will be taken up in 
this order. 


aCx"*") Equation 


The modification to the equations involved in the development of the 
A(x"^) equation are listed as follows: Add to the right-hand side of the 

equation listed the term following. 


(A 3 T)j first equation Fy^ 
(A 3 T)^ second equation 
(A 39 )^ first equation Cy^ 

(A 39 )j second equation 


(84) foei'P 


where fr 


Cpg + iCp.3 


(91) [iA - 


If vx"*" is small, the added term in equation (91) is approximated by 

(51) (1,1-1,,..) 

Hence, if the missile does not roll (v = O), the force due to asymmetry 
will cause the transverse displacement. A, to increase with the square of 
the distance, x”^. 


|(x+) Equation 


The modifications to the equations involved in the development of 
the ^(x"*") equation are listed as follows: Add to the right-hand side 
of the equation listed the term following. 


^The h(x'‘') motion is affected also but the changes required to the 
il(x'^) equations will be omitted from this treatment since the |(x'^) 
motion is the one ordinarily reduced in the analysis of flight data. 
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(a4o) 

PSCyo 

2 

(a4i) 

pSCzo 

2 

(a43) 

pSICuiQ 

2 

(a44) 

pSZCno 

2 

(TT) 

Gy 

iO 

2 

(78) 

"2b 

2 

(80) 

2Ky^ 

( 81 ) 

P"no 

(85) 

fo6^^ 


wbere Tq 


Cp3 + iCp3 
2 


(87) moeiq> 

K^M3 " 

where mo = — 5 


The aerodynamic asymmetries will change the values of the coefficients 
d;j_ and dg of the integrated |(x+) equation ( 93 ) J the following term 
should be added to P(s), defined by equation (llO): 

S(Sfo - m^fo - z^mo) J 

where ^ is the variable of Integration replacing in functions of 

x"*" appearing in the integrand and 9 is considered to be a function of 
^ . The added term may be readily derived in the integration of the 



MCA TN 3350 


95 


simultaneous differential equations (85) and (87) if the integration is 
carried out by the Laplace transform and use is made of the convolution 
of the "natural” frequency and the "asymmetrical forcing function." A 
clear explanation of the operational mathematics involved is given in 
chapters I and II of reference (16) (particularly section l 4 on the 
convolution; Churchill states that the convolution is also known as the 
Faltung integral) . 

It may be of interest to note that if the roll rate vith respect to 
distance is constant, that is, if 

cp = VX+ + q>Q 

then the term added to P(s) becomes 

S(Sfp - mT^fp - 2^mo)e^^o 

z. (iv - S) 

I 

Hence if v = S^i or then (iv - S) = or ordinarily a small 

value. Consequently, the value of the corresponding term will become 
very large. In other words, the oscillation experiences a divergent 
resonance as the roll rate approaches the pitch or yaw rate.^ 


3(iV-S)x+ _ J 


^This result has been noted elsewhere in the literature; for example, 
see reference 4 . 
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Figure 1.- Continued. 
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(c) Orientation of velocity vector with respect to stability axes. 


Figure 1 .- Continued. 



Figure 2.- Space geometry for projectiles and missiles with 90 
rotational and mirror symmetry. Orientation of body axes and 
pseudo-stability axes with respect to earth axes. 
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